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Abstract. We present a compilation of spectroscopic observations of the sgB[e] star CI Cam, the optical counter- 
part of XTE J0421+560. This includes data from before, during, and after its 1998 outburst, with quantitative 
results spanning 37 years. The object shows a rich emission line spectrum originating from circumstellar material, 
rendering it difficult to determine the nature of either star involved or the cause of the outburst. We collate all 
available pre-outburst data to determine the state of the system before this occurred and provide a baseline for 
comparison with outburst and post-outburst data. During the outburst all lines become stronger, and hydrogen 
and helium lines become significantly broader and asymmetric. After the outburst, spectral changes persist for 
at least three years, with Fell and [Nil] lines still a factor of ~ 2 above the pre-outburst level and He I, Hell, 
and Nil lines suppressed by a factor of 2-10. We find that the spectral properties of CI Cam are similar to other 
sgB[e] stars and therefore suggest that the geometry of the circumstellar material is similar to that proposed for 
the other objects: a two component outflow, with a fast, hot, rarefied polar wind indistinguishable from that of a 
normal supergiant and a dense, cooler equatorial outflow with a much lower velocity. Based on a comparison of 
the properties of CI Cam with the other sgB[e] stars we suggest that CI Cam is among the hotter members of the 
class and is viewed nearly pole-on. The nature of the compact object and the mechanism for the outburst remain 
uncertain, although it is likely that the compact object is a black hole or neutron star, and that the outburst was 
precipitated by its passage through the equatorial material. We suggest that this prompted a burst of supercritical 
accretion resulting in ejection of much of the material, which was later seen as an expanding radio remnant. The 
enhanced outburst emission most likely originated either directly from this supercritical accretion, or from the 
interaction of the expanding remnant with the equatorial material, or from a combination of both mechanisms. 

Key words. stars:individual(CI Cam) - starsiemission line, B[e] - stars:radio emission stars:CI Cam - binaries: 
close 



1. Introduction 

Send offprint requests to: R. I. Hynes, e-mail: 0n 1998 A P ril 2 Smith et aL ( l 1998 D reported an 
rih@astro.soton.ac.uk RXTE All-Sky Monitor (ASM) detection of a bright, 

Guest Observer, McDonald Observatory, University of rapidly rising X-ray transient desi gnate d XTE J0421+560 



Texas at Austin. Subsequently, Marshall et al. (1998) used the RXTE 
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Proportional Counting Array (PCA) to refine the best fit 
position with an error circle of I arcmin radius. The bright 
(V ~ II) B[c] star CI Cam(=MWC 84) was found to lie 
near to the centre of this error circle. Spectroscopic obser- 
vations by Wagner et al. ( |1998| ) on 1998 April 3 revealed 
a rich emission line spectrum, similar to that reported by 
Downes (1984; see Sect. |^), but with the presence of Hen 
emission features. These features had not been reported 
in previous spectra, and so by analogy to other X-ray bi- 
naries Wagner et al. (1998) proposed it to be the optical 
counterpart of XTE J0421+560. Photometric observations 



of the sourc e at this time (e.g. R obinson et al. [L998| , Garcia 
et al. 1998 , Hynes et al. 1998 ) showed that CI Cam was 
some 2-3 mag brighter than ha d previ ously been reported. 

Hjcllming & Mioduszewski ( 1998a ) reported the detec- 
tion of a transient 19mJy radio source at 1.4 GHz, cor- 
responding to the optical position of CI Cam on 1998 
April 1, thus confirming the identification of CI Cam 
as the optical counterpart. Subsequent observation of 
rapid radio variability established that the radio emis- 
sion was of non-thermal (synchrotron) origin (Hjellming 
& Mioduszewski 1998b). Long term observations indicate 
that after the initial flare the radio emission underwent 
an unusually slow decay, with a 15 GHz flux of ~ 1.5 mJy 
about 40 months after the initial outburst (Pooley, priv. 
comm.). High spatial resolution maps obtained after the 
outburst indicated the presence of a clumpy ejection neb- 
ula (Mioduszewski et al., in preparation). These ejecta 
expand at ~ 1.0 — 1.5 mas d -1 , corresponding to an ex- 
pansion velocity ~ 5000 km s -1 for a distance of 5kpc. 

Given the distance estimates (and optical luminosity 
impli ed) for CI Cam (e.g. log L/Lg> 4.86; Clark et al. 
2000| ; Robinson, Ivans & Welsh |2002j) it is clear that if 
it is a binary, as is likely, it is a high mass X-ray bi- 
nary (HMXB). However CI Cam does not sit comfortably 
within the traditional divisions of HMXB mass donors 
into classical Be stars (~ 70 per cent) and OB supergiants 
(~ 30 per cent). While its luminosity suggests that it be- 
longs to the later subset, such systems are typically short 
period binaries which accrete via Roche lobe overflow or 
direct wind fed accretion producing persistent X-ray emis- 
sion (typically modulated at the orbital period) . The pres- 
ence of a rich emission line spectrum including forbidden 
lines, and a near IR excess due to hot dust (the observa- 
tional criteria for the B[c] phenomenon; Allen & Swings 



1976; Lamers et al. 1998) also mark a distinction from the 



other supcrgiant HMXB systems. Among the stars show- 
ing the B[e] phenomenon, the high luminosity of CI Cam 
makes it a Galactic counterpart to the Magellanic Cloud 
supergiant B[e] stars (sgB[e] stars) and we will refer to it as 
such for the rest of this work. CI Cam therefore appears 
to be the first bona fide sgB[e] star HMXB known, al- 
though direct evidence for binarity has proven elusive and 
a chance encounter of a compact object with the sgB[e] 
star although highly unlikely, cannot be ruled out. 

This work presents a compilation of spectroscopy ob- 
tained before, during and after the 1998 outburst. A com- 
panion photometric compilation has been presented by 



Clark et al. (200C). Some of the outburst data included 
here has previously been presented by Barsukova et al. 
(|1998| ) and Barsukova et al. Q2002j ). In Sect. | we sum- 
marise the available pre-outburst data, including archival 
data with quantitative spectroscopy spanning ~ 30 years 
before the X-ray outburst and additional unpublished pre- 
outburst spectra. From these we identify the typical pre- 
outburst strengths of spectral lines and discuss their sta- 
bility. We then in Sect. || describe a series of new spec- 
tra running from a few days after the X-ray outburst to 
~ 3 years later. In Sect. |we discuss extinction and dis- 
tance estimates for the system and in Sect. |5 we review 
what is known about the mass donor star. Sect. 6J examines 
the changes in the continuum flux distribution and Sect. [?] 
the spectral lines and how these evolve through the out- 
burst. Sect. [| tests for the presence of shorter timescale 
variability. Finally in Sect. ^ we will discuss how all of 
these clues can help us build a picture of the nature of the 
system and the outburst mechanism and in Sect, [l^ we 
summarise our conclusions. 



2. Pre-outburst observations 

Since CI Cam is bright and has a strong emission line 
spectrum, there exists a relatively large data set dating 
from 1933 to the present day. We present a summary of 
these observations, which include previously unpublished 
optical spectra from 1987 and 1994, to quantify an average 
'quiescent' state for CI Cam. Table || summarises the pre- 
outburst equivalent widths. 

2.1. Published archival spectroscopy 

The first obser vations of CI Cam were made by Merrill & 
Burwcll ( 1933| ) and showed strong Hi Balmer and weak 
Fe II emission. He I was seen to be strongly in emission 
but Hell was absent. The first quantitative spectr oscopi c 
observations of CI Cam date from Chkhikvadze ( 1970 ). 
Spectra obtained during the period 1964-67 show strong 
H I and He I emission lines with no evidence of He n emis- 
sion. Fen lines were also in emission, the strongest being 
multiplets 27, 28, 37 and 38; there was no evidence for the 
Mgn 4481 A line. Based on the s trength of the Balmer 
emission lines Chkhikvadze ( 1970] ) suggested a spectral 
class of 06-BO for C I Cam . 

Allen & Swings ( 1976| ) describe further spectroscopic 
observations of CI Cam (although no date for the obser- 
vations is provided). They note numerous emission lines 
of He I, Fen and Sin. [Nil] lines at 5755 and 6584 A were 
present (the latter in the wings of the very strong Ha line) ; 
[Oi] lines at 6300 & 6363 A are possibly also detected. An 
estimate for the density of the circumstellar envelope of 
N c > 10 6 — 10 7 cm~ 3 was derived from the strength of 
the [Nil] lines. Higher excitation lines such as [Om] and 
[Nm] appeared to be absent. 

A 4000-7000A spectrum was obtained by Downes 



(|1984J) in 1984 January and once again was dominated 
by strong Hi, Hei and Fen lines. 
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Miroshnichcnko (1995) describes observations made 
between 1986 September and 1987 December. The spec- 
trum was dominated by strong H I and He I lines, with nu- 
merous weak Fe n lines also present (as was C n emission 
at 4267 & 7234A). The Hi Balmer lines were symmetrical 
and single peaked, with wings extending to ±250 km s -1 
for Ha and H/3. 

Finally, Jaschek & Andrillat ( [2000 ) report observa- 
tions from 1992 and 1998, the latter only two months 
before the outburst. 450 emission lines were observed, 
~ 55 percent from Fen, the remainder including Hi, Hei, 
Oi, Ni, Sin, Mgn, [Oi], [Nil], and [Fen]. The Balmer 
lines show a very steep Balmer decrement. Oi 8446 A is 
unusually strong, likely due to hy/3 fluorescence. Very few 
differences between the 1992 and 1998 epochs were seen. 
Since the latter spectra represent the highest quality pre- 
outburst spectra, as well as the closest to the X-ray out- 
burst, we reanalysed them to establish a quantitative pre- 
outburst baseline. 



2.2. Unpublished spectra taken at the Special 
Astrophysical Observatory 

Three further unpublished observations of CI Cam were 
obtained by A. Miroshnichenko at the 6 m BTA telescope 
of the Special Astrophysical Observatory of the Russian 
Academy of Sciences (SAO RAS) on 1987 April 5 (-4460- 
5000A) and 1994 January 21 (-4000-4950 & -5700- 



The observations were made with a medium-resolution 
spectrograph SP-124 and a 1024-elemcnt one-dimensional 
photo-electric detector (Drabek ct al. 1986 ). Data reduc- 
tion was performed with the SIPRAN software developed 
at SAO RAS; no attempt at flux calibration was made. 
Due to the poor dynamic range of the scanner it was pos- 
sible to obtain either profiles from the strong lines on an 
underexposed continuum, or a properly exposed contin- 
uum with saturated strong lines. Of the data presented 
here, the 1987 spectrum was exposed to search for weak 
lines, while the two spectra obtained in 1994 were ob- 
tained to determine the line profiles of the stronger emis- 
sion lines. Consequently the continuum (and weaker lines) 
were poorly exposed for these spectra. We have chosen to 
include an equivalent width for the lines identified from 
the 5700-6700A spectrum, although we caution that the 
continuum level was uncertain in this spectrum. 



2.3. Archival spectra from Haute-Provence 

We obtained the 1998 spectra reported by Jaschek & 



Andrillat (2000) from the archive of the Observatoire de 
Haute-Provence (OHP). These spectra were obtained with 
the Aurelie spectrograph on the 1.52 m telescope using a 
Thomson array of 2026 photodiodes and a 300 line mm -1 
grating, yielding a resolution of about 1.3 A. 
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Table 1. Pre-outburst equivalent widths, in A, of se- 
lected strong emission lines. Minus signs have been 
dropped for convenience as all lines are in emission. Hei 
4921 and 5015 A have been omitted as these are unre- 
solved from strong Fen l ines. Sources ar e 1 C hkhikvadze 
( |1970| ), 2 Mirosh nichenko (|l995|) 3 Section |J, 4 Jaschek & 
Andrillat ( [2000| ) and Section 2JJ. Lines with no value given 
are not reported, outside the spectral range or too weak to 
reliably measure. Bracketed measurements were obtained 
from spectra with a poorly exposed continuum rendering 
the equivalent widths uncertain. 



2.4. The 'average' quiescent spectrum 

We plot the pre-outburst optical spectrum (from 1998 
OHP data) in Fig. %i and in Table [l] we compile the equiv- 
alent widths of the strongest emission lines in quiescence. 
Coverage is clearly patchy, but no line exhibits dramatic 
variations in equivalent width. In particular the relative 
strengths of Balmer and He lines remain roughly the same. 
Downes ( 1984 ) does not include equivalent widths, but the 
relative line strengths plotted in his spectra (e.g. H7 to 
Hei 4471) look very similar to the results in Table [j], and 
also to Fig. |l|a. It therefore appears that the pre-outburst 
spectrum was stable for at least — 30 years before the 
outburst, although large gaps in coverage are obviously 
present. As discussed in Sect. 0, the post-outburst spec- 
trum is clearly different to this, but it remains to be seen 
whether the differences represent a delayed return to the 
pre-outburst quiescent state, or a new different quiescent 
state. 



3. Outburst and post-outburst observations 

We next present a series of spectroscopic observations ob- 
tained during 1998 April and May, covering the outburst 
and immediate aftermath. Spectroscopy taken during this 
time indicates the presence of high excitation lines of He 11 
at the time of the outburst which fade on a timescale 
— 10 days; similar variability is seen in most other lines 
observed during this period. Such behaviour occurs over 
longer timescales than the X-ray outburst; XTE observa- 
tions reveal that the X-ray flux had already d roppe d by a 
factor of a hundred by April 4 (Belloni et al. 1999 ). 
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Fig. 1. Comparison between normalised spectra from a) before the outburst (1998 January 27-28), b) early in the 
outburst (1998 April 4), and c) after the outburst (2000 February 6 for A > 7100 A and 2000 December 1 elsewhere). 
Major features are labelled. The strongest groups of Fell lines are also indicated, but other species are present in these 
complexes and Fe n lines are also found elsewhere. 



We obtained both blue and red spectra during the out- 
burst, thus enabling us to follow the evolution of the out- 
burst and determine changes in the circumstellar environ- 
ment during this period. Spectra were obtained from four 
telescopes in the month following the X-ray outburst: the 
1.5 m at the FLWO, the 6 m telescope of the SAO RAS, the 
4.2 m WHT on La Palma and the McDonald Observatory 
2.7 m. Subsequent observations obtained during the pe- 
riod between 1998 May and 2001 October constrain the 
long term post-outburst behaviour of the system. All of 
these observations are summarised in Table [2] and more 
details are given in the following subsections. 



3.1. F. L. Whipple Observatory observations 

A series of spectra were obtained using the 1.5 m tele- 
scope at the F. L. Whipple Observatory (FLWO) and the 
FAST spectrometer. These observations began soon af- 
ter the outburst was discovered and continued until the 
end of 1998. The spectra were obtained with two different 
gratings, of 300 and 1200 lines mm -1 , giving 3.0 and l.lA 
resolution respectively. 

These observations included a series of 370 2 s spectra 
with ~ 1.1 A resolution and a dispersion of 0.4A pix , 
covering 6040-7035 A which were obtained on the night of 
1998 April 2. We searched for shifts in the velocities of the 
emission lines by cross-correlating all the spectra against 
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Date Telescope UT Start UT End Number Wavelength Resolution 

Range(A) (A) 
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Table 2. Summary of spectroscopic observations, both pre-, during, and post-outburst which were quantitatively 
analysed for this work. In some cases multiple spectra from the same night covering the same wavelength range have 
been merged on one line; in this case start and end times refer to the first and last spectrum and do not indicate total 
exposure time. Continued on the next page. 



a single template. The largest shifts are ±12kms -1 , the 
rms in the velocities is 3.6 km s . Given the high airmass 
(from 1.5 to 2.0) and the 2.0arcsec slit, these variations 
are likely due to variable illumination of the slit rather 
than intrinsic motions in the emission lines, as the velocity 
resolution was about 50kms _1 . Nonetheless, we searched 
for periodicities in the velocities using the algorithm of 
Stcllingwerf ( 1978 ), but found nothing significant. We will 
further discuss evidence for variations in line strengths and 
profiles in Sect. 



Photometrix 1024x1024 CCD detector and the Bl grat- 
ing (600 lines mm -1 ) providing a spectral dispersion of 
2.4Apixel~ 1 in two overlapping spectral ranges 3700- 
6100 and 5000-7500 A. The slit width employed was either 
1 arcsec or 2 arcsec, depending on seeing conditions, pro- 
viding a spectral resolution of ~ 4 A or ~ 8 A respectively. 
Primary reduction of the CCD spectra including correc- 
tions for the dark current, debiasing and sky subtraction 
was accomplished with MIDAS (1996 Nov. version). He- 
Ne-Ar lamps were used for wavelength calibration. 



3.2. Special Astrophysical Observatory observations 

Outburst observations spanning 1998 April to May 
were made using the 6 m BTA telescope of the 
SAO RAS. A medium resolution spectrograph, SP- 
124 (at the Nasmyth-1 focus), was used with a 



3.3. WHT/UES observations 

CI Cam was observed at high resolution using the Utrecht 
Echelle Spectrograph (UES) on the William Herschel 
Telescope (WHT) on La Palma on 1998 April 9-11. 
Spectra were extracted using a simple (non-optimal) ex- 
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16/10/98 


FLWO 1.5-m 


10:44 


10:47 


2 


3650-7520 


3.0 


23/10/98 


FLWO 1.5-m 


10:49 


10:51 


3 


3580-7450 


3.0 


29/10/98 


FLWO 1.5-m 


09:03 


09:07 


2 


3650-7520 


3.0 


01/11/98 


Asiago 1.82-m 


02:20 


02:35 


1 


6360-8060 


3.3 


12/11/98 


FLWO 1.5-m 


10:02 


10:13 


3 


3630-7500 


3.0 


13/11/98 


FLWO 1.5-m 


11:41 


11:44 


3 


3630-7500 


3.0 


14/11/98 


FLWO 1.5-m 


10:02 


10:06 


3 


3630-7500 


3.0 


15/11/98 


FLWO 1.5-m 


08:14 


08:17 


3 


3630-7500 


3.0 


1 ez /1 1 /no 
16/11/90 


1LWU l.o-m 


1 t\.A 1 

1U:41 


1 (\.A A 

1U:44 


Q 
■i 


Q£?QA TKAA 

dodU-l oUU 


■3.U 


17/11/98 


FLWO 1.5-m 


11:08 


11:10 


4 


3630-7500 


3.0 


18/11/98 


FLWO 1.5-m 


09:06 


09:18 


4 


3630-7500 


3.0 


19/11/98 


FLWO 1.5-m 


09:40 


09:44 


3 


3660-7530 


3.0 


21/11/98 


FLWO 1.5-m 


09:31 


09:34 


3 


3660-7530 


3.0 


22/11/98 


FLWO 1.5-m 


08:35 


08:39 


3 


3660-7530 


3.0 


23/11/98 


FLWO 1.5-m 


09:04 


09:07 


3 


3660-7530 


3.0 


24/11/98 


FLWO 1.5-m 


10:24 


10:27 


3 


3660-7530 


3.0 


25/11/98 


FLWO 1.5-m 


07:56 


08:01 


3 


3660-7530 


3.0 


26/11/98 


FLWO 1.5-m 


07:47 


07:51 


3 


3660-7530 


3.0 


27/11/98 


FLWO 1.5-m 


08:23 


08:27 


3 


3660-7530 


3.0 


30/11/98 


FLWO 1.5-m 


11:17 


11:22 


3 


3660-7530 


3.0 



Table 2. (continued) 



traction with sky subtraction using echomop software. 
Cosmic rays were removed by hand. 

3.4. McDonald Observatory observations 

Intermediate dispersion spectroscopy was obtained on the 
nights of 1998 April 18-20 using the Large Cassegrain 
Spectrometer (LCS) on the McDonald Observatory 2.7m 
telescope. Each night involved a series of observations 
spanning about 2 hours to search for variability; this will 
be discussed further in Sect. |8|. Reduction and extraction 
were performed using standard iraf tasks. Observations 
of neon arc lines were performed at 20 min intervals to 
provide wavelength calibration. Corrections were made to 
the interpolated wavelength solution using night sky emis- 
sion lines. As for the FLWO observations, no significant 
wavelength variations were found. 



3.5. WHT /ISIS observations 

Intermediate and high dispersion spectroscopy was ob- 
tained on several nights from 1998-2001 using the ISIS 
dual-beam spectrograph on the WHT. Most exposures 
used the blue arm with the EEV10 or EEV12 CCD. 
The second spectrum from 2001 April 28 used the red 
arm with the TEK4 CCD. Slit widths were 0.7-1.0 arcsec. 
Reduction and extraction were performed using standard 
iraf tasks. Copper-neon and/or copper-argon lamps were 
observed to provide wavelength calibration. No flux cali- 
bration of these spectra was attempted. 

3.6. Loiano & Asiago Observatory observations 

Post-outburst spectra of the source were obtained in 1999 
January 5 and 2000 February 6 using the 1.52 m G. D. 
Cassini telescope at the Loiano Observatory (Italy). The 
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Date 


Telescope 


UT Start 


UT End 


Number 


Wavelength 
Range (A) 


Resolution 

(A) 


10/12/98 


FLWO 1.5-m 


05:59 


06:03 


3 


3660-7530 


3.0 


12/12/98 


FLWO 1.5-m 


07:29 


07:32 


3 


3660-7530 


3.0 


13/12/98 


FLWO 1.5-m 


05:29 


05:33 


3 


3660-7530 


3.0 


14/12/98 


FLWO 1.5-m 


09:25 


09:29 


3 


3660-7530 


3.0 


19/12/98 


FLWO 1.5-m 


06:07 


06:11 


3 


3660-7530 


3.0 


20/12/98 


FLWO 1.5-m 


06:22 


06:27 


3 


3660-7530 


3.0 


21/12/98 


FLWO 1.5-m 


07:35 


07:39 


3 


3660-7530 


3.0 


22/12/98 


FLWO 1.5-m 


06:08 


06:12 


3 


3660-7530 


3.0 


23/12/98 


FLWO 1.5-m 


06:11 


06:15 


3 


3660-7530 


3.0 


25/12/98 


FLWO 1.5-m 


08:32 


08:36 


3 


3660-7530 


3.0 


26/12/98 


FLWO 1.5-m 


04:52 


04:56 


3 


3660-7530 


3.0 


27/12/98 


FLWO 1.5-m 


05:14 


05:17 


3 


3660-7530 


3.0 


28/12/98 


FLWO 1.5-m 


03:14 


03:18 


3 


3660-7530 


3.0 


03/01/99 


OHP 1.52-m 


00:26 


01:55 


1 


8040-8900 


1.3 


03/01/99 


OHP 1.52-m 


03:30 


05:30 


1 


4060-4930 


1.3 


03/01/99 


OHP 1.52-m 


20:05 


21:35 


1 


4860-5730 


1.3 


03/01/99 


OHP 1.52-m 


22:11 


23:41 


1 


5560-6430 


1.3 


04/01/99 


OHP 1.52-m 


00:04 


01:34 


1 


6250-7110 


1.3 


04/01/99 


OHP 1.52-m 


02:33 


04:33 


1 


7060-7930 


1.3 


05/01/99 


Loiano 1.52-m 


21:17 


21:27 


1 


6360-8220 


3.3 


29/01/00 


Calar Alto 3.5-m 


23:44 


23:49 


1 


3700-6800 


6.1 


06/02/00 


Loiano 1.52-m 


21:11 


21:41 


1 


3500-5400 


5.5 


06/02/00 


Loiano 1.52-m 


20:33 


21:06 


2 


3530-8830 


8.3 


19/07/00 


Skinakas 1.3-m 


02:30 


02:35 


1 


5550-7550 


4.5 


16/10/00 


Skinakas 1.3-m 


01:37 


01:40 


1 


5250-7300 


4.5 


01/12/00 


WHT 4.2-m 


20:22 


20:37 


3 


3700-7200 


3.6 


01/12/00 


WHT 4.2-m 


21:06 


21:37 


3 


3870-4310 


0.3 


28/04/01 


WHT 4.2-m 


20:53 


21:03 


1 


3600-4030 


0.3 


28/04/01 


WHT 4.2-m 


20:45 


21:03 


15 


6300-6700 


0.6 


08/08/01 


Skinakas 1.3-m 


01:42 


01:47 


1 


5460-7430 


4.5 


13/09/01 


Skinakas 1.3-m 


00:45 


00:48 


1 


5230-7200 


4.5 


23/10/01 


Bok 2.3-m 


10:30 


10:46 


3 


3885-5030 


1.8 



Table 2. (continued) 



telescope was equipped with the Bologna Faint Object 
Spectrograph and Camera (BFOSC) and a 2048x2048 
Loral CCD. In January 1999 grism #8 was used, giv- 
ing a resolution of ~ 3A, and in February 2000 the low- 
dispersion grisms #3 (blue) and #4 (red) were used. One 
further spectrum was taken on 1998 November 1st us- 
ing the 1.82m telescope at Asiago Observatory (Italy), 
equipped with the Asiago Faint Object Spectrograph and 
Camera (AFOSC) and a 1024x1024 thinned SITc CCD. 
This instrument is identical to BFOSC and grism #8 was 
used, giving the same resolution. 



3. 7. Haute Provence observations 



3.8. Calar Alto observations 

One low-resolution flux-calibrated spectrum was taken on 
2000 January 29 using the blue arm of the TWIN spec- 
trograph on the 3.5 m telescope at Calar Alto (Spain) and 
a 2000x800 thinned SITe CCD. 

3.9. Skinakas Observatory observations 

Spectra taken on 2001 July 19 & October 16 were 
taken with the 1.3 m f/7.7 Ritchey-Chretien telescope at 
Skinakas Observatory (Crete, Greece). This was equipped 
with a 2000x800 ISA SITe CCD and a 1302 lines mm" 1 
grating (using a 80/im slit), giving a dispersion of 
—1 Apixcl -1 . 



Post-outburst spectra obtained from OHP in 1999 used 
the same configuration as the 1998 pre-outburst obser- 
vations, hence the details are the same as described in 
Sect.Es. 



3.10. Steward Observatory observations 

Finally three spectra were obtained using the Bok 2.3 m 
telescope of the Steward Observatory with the Boiler and 
Chivens long slit spectrograph and a 1200x800 CCD. A 
first order 1200 lines mm -1 grating was used and there 
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were no filters. Flux calibration was not reliable as there 
were significant slit losses. 

4. Interstellar extinction and distance estimates 



There remain large uncertainties about the interstellar ex- 
tinction to CI Cam, and even more so concerning the dis- 
tance. A number of extinction estimates have been made 
but these are complicated by the uncertainty about how 
much of the measured extinction is intrinsic to the source; 
different methods may be more or less sensitive to this de- 
pending on what is actually measured. Consequently it is 
of value to compare as many independent methods as pos- 
sible. We discuss only recent estimates as these are likely 
to use higher quality data and/or more sophisticate d mod - 
els, superseding earlier estimates (e.g. Chkhikvadze 1970|) . 

A number of attempts have been made to fit the broad 
band optical-far IR spectral energy distribution (SED) 
using a model including extinction as a free parameter. 
Belloni et al. ( 1999 ) fitted the SED using a Kurucz at- 
mosphere and a dust model and derived E(B — V) = 
1.18 ± 0.04mag and hence Ay ~ 4.4mag. Clark et al. 
(200C) used a similar model to derive Ay ~ 3.71 mag. 



Zorcc (199S) used a model in which the interstellar ex- 
tinction and distance were constrained to be consistent 
rather than independent, and estimated d ~ 1.75 kpc, 
E(B - V)i sjn ~ 0.8 mag and ^4y,i cai ~ 2.4 mag, imply- 
ing a total extinction somewhat higher than the other two 
estimates. 

Orlandini et al. (fOOFj) used emission line ratios to esti- 
mate E(B-V) = 1.54 mag (from He lines) or E(B — V) = 
1.02 mag (from H lines). This is based on theoretical pre- 
dictions of the ratios of the strongest lines in CI Cam. In 
Be stars and other early OB type stars with dense winds, 
and likely also in sgB[e] stars, these lines are subject to 
effects of non-local thermodynamic equilibrium (NLTE) 
so this method will not be very reliable. 



Robinson et al. (2002) use the 2175 A interstellar ab- 
sorption feature in post-outburst HST data to measure 
E(B -V) = 0.85 ± 0.05 mag and A v = 2.3 ± 0.3 mag; the 
large error on Ay reflects the uncertainty in the choice of 
extinction curve parameterised by Ry = Ay /E(B — V). 
We can attempt to estimate the Na D equivalent width, 



although as noted by Munari & Zwitter (1997) this is 
insensitive for reddened objects as it tends to saturate. 
It is also complicated by the presence of Na D emission. 
Fortunately, our echelle spectra obtained in outburst re- 
solve the NaD components as sharp absorptions within 
the emission line (Fig. ||a). Reconstructing the unabsorbed 
line profile is obviously uncertain, but assuming it is sym- 
metric we can estimate that the equivalent widths of the 
two NaD2 (5890 A[]) components observed are 0.36 and 
0.52 A. The latter component does appear saturated, so 
this only provides a lower-limit to the reddening. If the 



1997 



incorrectly label the 



1 We note that Munari & Zwitter 
5890A line as NaDl, and 5896A as NaD2. Their work is, 
however, internally consistent; Munari, priv. comm. 




5885 5890 5895 

Wavelength (A) 



5900 




-60 -40 -20 

LSR velocity (kms -1 ) 



20 



Fig. 2. a) NaD line profiles observed during outburst 
(1998 April 9). The two interstellar absorption compo- 
nents are clearly visible; the longer wavelength one is sat- 
urated. The dashed line shows the reconstructed NaD2 
line profile assuming it is symmetrical, b) Closeup of the 
absorption components indicating the velocities with re- 
spect to the LSR. The structure in the two lines is con- 
sistent. In this direction the LSR velocity is expected to 
become more negative with distance, so the distance in- 
creases to the left. A quantitative distance scale depends 
on the assumed rotation curve of the Galaxy. 



profile is actually asymmetric with the same extended 
blue wing seen in many other lines then this will in- 
crease the inferred equivalent width, implying a redden- 
ing further above our lower limit. Using the calibration 
of Munari & Zwitter ( 1997| ) for each of the two compo- 
nents we derive a lower limit for the combined reddening 
of E(B-V) £0.5 mag. 

Our spectra also reveal several diffuse interstellar 
bands (DIBs). Our principal reddening cali brato r is the 
A5780 band, using the calibration of Herbig ( 1993j ). Using 
the WHT /UES spectra we measure an EW of 320 ± 30 mA 
implying E(B - V) = 0.64 ± 0.06 mag. The DIB at 
5797 A is too blended with emission lines to reliably mea- 
sure but we can approximately estimate 140 ± 30 mA im- 
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plying E(B - V) = 0.9 ± 0.2 mag (Herbig |1993| ). For 
5705 A we measure 110 ± 20 m A implying E(B - V) = 
0.49 ± 0.09 mag (Herbig 1975 ). We also use an average 
of several WHT/ISIS and Skinakas spectra to estimate 
an equivalent width of 180 ± 30 mA for A6203 implying 
E(B - V) = 0.62 ± 0.11 mag (Herbig [1975D , and a cen- 
tral depth of 0.90-0.97 for the very broad A4428 featur e 
implying 0.3 < E(B - V) < 1.1 (Krelowski et al. |1987p 



The 6284 A line may also be present but is complicated by 
Telluric absorption. All of the reliable DIB measurements 
are consistent with E(B — V) ~ 0.65 mag as suggested by 
our principal calibrator, 5780 A. 

To summarise, methods based on interstellar absorp- 
tion features (2175 A, NaD, DIBs) all suggest E(B-V) ~ 
0.5 — l.Omag, and hence visual extinction Ay ~ 1.5 — 
3 mag. Methods based on fitting the SED or line ratios, 
however, favour larger total extinction Ay ~ 4 mag. This 
could arise if the absorption features are only produced 
in the interstellar medium, and not by local dust extinc- 
tion, as a consequence of different chemical compositions 
and/or grain sizes. The SED methods are all model de- 
pendent, however, and it may be that the models used are 
either incorrect or incomplete. This leaves us with some 
ambiguity about how to deredden spectra of CI Cam. The 
higher values inferred from SED modelling are most rele- 
vant, in the sense that they directly measure the distortion 
of the broad band spectrum, but they are also model de- 
pendent. Measurements based on absorption features are 
more objective, but do not measure the same thing. 

Distance estimates for CI Cam are plagued by even 
more uncertainty than extinction estimates. A number 
of wo rks ha ve der ived distances of 1- 2 kpc (Chkhi kvadz e 
1970| ; Zorec |l998[ Belloni et al. |l999| ; Clark et al. |2000|) . 



All of these, however, involve uncertain assumptions about 
the luminosity class of the star, or the relation betwe en in - 
terstellar absorption and distance. Robinson et al. ( 2002 ) 
challenged this conclusion, arguing for a larger distance. 
Based on spectroscopic similarities to the largest and most 
luminous sgB[e] stars they concluded that the distance 
had to be larger than 2 kpc. Since the line of sight passes 
well above the warped Galactic plane for 2-6 kpc, they 
argued that a young object like CI Cam must lie beyond 
rather than within this distance range. In support of this 
they note that its radial velocity is consistent with a dis- 
tance of ~ 7 kpc, assuming differential Galactic rotation. 
These arguments, however, make the assumptions that 
CI Cam is comparable to the most luminous sgB[c] stars, 
and that it does not have significant peculiar velocity. 

More objective constraints on the distance are harder 
to obtain. We can, however, make one estimate which 
is almost completely independent of what CI Cam is. 
As already described, the NaD lines show sharp ab- 
sorption components. The velocities of these in the lo- 
cal standard of rest (LSR) are —4.5 ± 0.6 km s -1 and 
— 34.8±0.4kms~ 1 for the stronger and weaker components 
respectively (see Fig. ||b). The stronger component is sat- 
urated so it could itself involve multiple components with 
LSR velocities up to — lOkrns -1 . The weaker component 



may also have an additional component in its blue wing at 
around — 45kms~ 1 ; both the NaDl and NaD2 lines sug- 
gest such a feature at this velocity. All of these components 
are somewhat redshifted with respect to CI Cam, which 
has an LSR velocity of — 51kms _1 (Robinson et al. 2002). 
It is therefore unlikely they are associated with circumstcl- 
lar material; much more likely is that they are interstel- 
lar gas. If so then we expect them to be moving in the 
plane in near circular orbits following Galactic rotation, 
and their velocities can be used to estimate their distances, 
and hence a lower limit on the distance of CI Cam. There 
are obviously uncertainties introduced by non-circular mo- 
tions and an imperfectly known rotation curve. We esti- 
mate distances using the rang e of r otation curves illus- 
trated by Oiling & Merrifield ( 1998| ), with the spread in 
values giving an estimate of the uncertainty in the mea- 
surement. The saturated component is clearly associated 
with the Local Arm, with a maximum velocity — lOkms" 1 
corresponding to a maximum distance ~ 1 kpc. We might 
expect the next feature to correspond to the Perseus Arm 
at a distance ~ 2.5 kpc, but the Perseus Arm is not well de- 
fined in the direction of CI Cam as our line of sight passes 
well above the Galactic plane at that point (c.f. discus- 
sion by Robinson et al. 2002), so it is unsurprising that it 
produces no absorption feature. The next feature out, at 
— 35kms _1 , is at an implied distance of 3.9-5.6kpc, and 
probably corresponds to the next spiral arm out, where 
our line of sight passes back into the warped outer disc; 



Cam OB3 (Humphreys et al. 1978) likely belongs to this 
arm. If the third weak feature, at — 45kms _1 , is real and 
also indicative of Galactic rotation, then it suggests an 
even greater minimum distance of 6-8 kpc, dependent on 
the assumed rotation curve. 

An alternative empirical appr oach is to use the radial 
velocity maps of Brand & Blitz ( 1993 ), derived from Hn 
regions and reflection nebulae. This will take account of 
non-circular motions. If anything, these suggest an even 
larger distance. The coverage in the direction of CI Cam is 
sparse, but representative distant objects from their sam- 
ple (S208, S211 and S212) have LSR velocities between 
-30 and -38kms _1 and distance estimates of 5.9-7.6kpc. 

The distance implied for CI Cam by the interstellar 
features is thus large; it is at least 4 kpc, and may well 
be beyond 6 kpc. This conclusion is essentially consistent 
with that derived independently by Robinson et al. (2002) 
using different methods. 

If the closest distance is adopted it is possible that 
CI Cam could be associated with Cam OB3, possibly as a 
runaway object. It lies ~ 2?7 from the centre of the associ- 



ation; the members identified by Humphreys (1978) lie at 
up to L°2 from the centre. The distance modulus adopted 



by Humphreys (1978) was 12.6 (3.3 kpc). CI Cam there- 
fore lies somewhat outside the association, and appears 
to be further away, but we cannot rule out the possibility 
that it is a runaway from Cam OB3. Alternatively it could 
lie beyond the spiral arm containing Cam OB3. 

Where it is necessary to assume a value for the redden- 
ing in what follows, we either take E(B — V) = 1.3 mag 
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(A v = 4 mag) or a range of 0.65 < E(B - V) < 1.4 
(2.0 < Ay < 4.4). For the distance, we follow Robinson et 
al. (2002) in adopting 5 kpc as a representative estimate, 
although we agree that it could be somewhat larger than 
this. Where a lower limit is more appropriate we assume 
d > 4 kpc. 

5. The nature of the components of CI Cam 

CI Cam clearly shows the observational characteristics of 
the Bjel phenomenon (Allen & Swings 1976; Lamers et 
al. 1998 ): strong Balmer emission lines (e.g. Fig. [j]), low 
excitation permitted lines (e.g. Fi g. |l|) , optical forbidden 
lines of [Fen] and [Oi] ( e.g. F ig. |12|a) and a strong in- 
frared excess (Clark et al. |2000D ; indeed CI Cam has long 
been considered a B[e] star a nd was included in the sam- 
ple of Allen & Swings |1976| . In view of the distance es- 
timates discussed above CI Cam clearly falls within the 
sgB[e] sub-class. The primary characteristic of ansgB[e] 
star, a luminosity of above 10 4 L Q (Lamers et al. 1998 ) is 
satisfied for any distance above 0.8 kpc and for the lower 
limit of 4 kpc that we have argued for above, this rises to 
10 5 4 L Q . Robinson et al. (2002) also note spectroscopic 
similarities to the most luminous sgB[e] stars, consistent 
with this. The other known types of stars showing the 
B[e] phenomenon are the pre-main sequence Herbig Ae/Be 
(HAeB[e]) stars, compact planetary nebula B[e] (cPNB[c]) 
stars and symbiotic B[e] (symB[e]) stars. HAeB[e] stars 
typically have luminosities of ^ 1O 45 L0, are associated 
with star forming regions and sometimes show evidence 
of infall, e.g. inverse P Cygni profiles. CI Cam is more lu- 
minous than this and exhibits a steep d ecrea se in the IR 
flux longward of ~ 10 /im (Clark et al. |2000|) , indicative 
of the absence of the cold dust that would be expected 
around an HAeB[e] star. Lamers et al. ( |l998|) did suggest 
that CI Cam could be a cPNB[e] star, but as these objects 
have lower luminosities (^ 10 4 L Q ) this is also ruled out 
by the distance estimates discussed in Section |[ SymB[e] 
stars are binaries also containing a cool giant which is usu- 
ally seen in the red or infrared spectrum. Some authors 
ha ve ide ntified CI Cam as a symbiotic (e.g. Barsukova et 
al. [2002 ), but late type features are not seen in CI Cam 
(except for the report of Miroshnichenko 1995 which was 
not corroborated by any other observations.) We believe 
that both the observed spectrum and the high luminosity 
therefore identify CI Cam as an sgB[e] star, as argued by 
Robinson et al. 



2002 



The sgB[e] class itself includes a range of objects with 
likely spectral types spanning B0-B9, and luminosities 
from 1O 4 2 -1O 61 L (Lamers et al. |1998[ ). It is therefore 
of interest to attempt a more precise classification for 
CI Cam. We might hope that the detection of photospheric 
lines from the sgB[e] star would be of great help. In the 
highest resolution post-outburst spectra of CI Cam, broad 
absorption wings are seen around higher order Balmer 
lines, most prominently H<5 and He (Fig. ^|). These may 
well be photospheric absorption lines from the sgB[e] star. 
Unfortunately we cannot obtain any detailed diagnostics 
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Fig. 3. Absorption wings surrounding H<5 and He lines, 
from 2000 December 1 WHT data. Unmarked lines are 
mostly Fell. 



from these lines because of the heavy contamination we see 
in the wings of the lines (both in terms of other lines and 
also probably residual wind emission); hence it is impos- 
sible to quantify the underlying photospheric spectrum. 
Indeed, the photospheric spectrum may be poorly defined 
for a star with a very high mass loss rate for which a bona 
fide photospheric radius is hard to determine. 

While the spectral type of the star is difficult to de- 
termine from the optical spectrum due to heavy contam- 
ination by the circumstellar material, the detection of 
P Cygni profiles in the UV resonance lines with absorp- 
tion to ~ 1000 km s -1 (Robinson et al. 2002 ) is more use- 
ful. By analogy with other sgB[e] stars we expect that the 
UV P Cygni profiles arise from a hot, polar wind which 
is similar to that of normal hot supergiants. The outflow 
velocity implied for CI Cam is directly comparable to that 
seen for other early B supergiants. The Si iv 1394, 1402A 
doublet functions as a powerful probe of temperature and 
luminosity for B stars (Walborn, Parker, & Nichols 1995). 
Early-B supergiants show a strong P Cygni wind profile, 
which evolves into pure absorption for mid-B stars and 
is absent in late B stars of all luminosity classes - the 
P Cygni profiles therefore suggest an early B classifica- 
tion. A similar conclusion can be drawn from the presence 
of P Cygni profiles for the Civ 1549,51 A doublet which 
is seen for supergiants of type B4 or hotter (and is seen 
in absorption for hot BO-2 dwarfs and giants and cooler 
supergiants). 

By comparison of the optical emission lines with 
other sgB[e] stars and luminous blue variables (c.f. 



Miroshnichenko 1996), we can further narrow the spectral 
type to B0-B2, as earlier spectral types show stronger 
He II emission than the very weak pre-outburst feature 
that we see, and later types show He I in absorption rather 
than emission. 

To summarise, the 'normal' star is an sgB[e] star, 
with likely spectral type B0-B2 and a luminosity of at 
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least 10 5 4 Lq, placing it among the hotter, more luminous 
sgB[e] stars (c.f. Robinson et al. 2002). 

The nature of the compact object implicated in the 
outburst is even more uncertain than that of the nor- 
mal star. It is widely assumed from the X-ray and 7-ray 
outburst observed that this must be a black hole or neu- 
tron star (e.g. Belloni et al. 1999; Robinson et al. 2002| ). 
Orlandini et al. (2000), however, have argued for a ther- 
monuclear runaway on the surface of a white dwarf. If 
the 4kpc lower limit on the luminosity is correct, how- 
ever, then the X-ray outburst was extremely luminous, 
Lx <; 2 x 10 38 ergs -1 . This corresponds to the Eddington 
limit for a neutron star, and for larger distances then 
a black hole becomes more likel y. Ho wever, neith er the 

Ueda et al 



X-ray spectral sh ape (O rr et al. [1998 



1998 



Revnivtsev et al. |l999j Belloni et al. 1999) no r the lack 
of rapid X-ray variability (Frontera et al. 1998 ; Belloni et 
al. 1999 ) are typical of accreting black holes or neutron 
stars. Equally, while there may be a hint of a flaring soft 
component (Frontera et al. 1998 ; Ueda et al. 1998|) , the 
X-ray spectrum is on the whole quite hard, and not domi- 
nated by a supersoft component as might be expected for 
a high luminosity accreting white dwarf. Consequently the 
identification must remain uncertain. 





-20.0 






'n 
X 


-20.5 - 


'(/) 


-21 .0 : 


cm 


-21 .5 : 


01 




CD 


-22.0 - 






Li_ 




cn 


-22.5 : 



_i 






-23.0 1 




14.60 14.70 14.80 

Log v 



14.90 



Fig. 4. Flux calibrated spectra obtained from FLWO early 
in outburst and after it. These have been dereddened 



using the Fitzpatrick (1999) extinction curve assuming 
Ay — 4.0. In outburst, the lines become much stronger 
and broader, Balmer jump emission appears, and the con- 
tinuum becomes redder. Note that these are long exposure 
spectra to ensure that the continuum is well defined, so 
the stronger emission lines are saturated. 



6. The spectral flux distribution 



Clark et al. (2000) considered the broad band spectral 
energy distribution from photometric measurements and 
examined colour changes during the outburst. They found 
that the source became systematically redder when it was 
brighter. Given the large contribution of emission lines to 
the spectrum, however, it is useful to attempt to isolate 
continuum changes, and quantify the line contribution. 

In Fig. |we show the first outburst spectrum obtained, 
which fortunately was flux calibrated, together with a 
comparable post-outburst one. It is clear that the red- 
der spectrum seen in outburst was not just due to line 
changes, but that the underlying continuum is redder. We 
will further discuss the origin of this enhanced continuum 
emission in Section 9.3. There is also an enhancement at 



higher frequencies as the Balmer jump appears in emis- 
sion, which explains why the (U — B) ind ex remained con- 
stant during outburst (Clark et al. 200(j| ) . 

In principle, we could also use flux calibrated spec- 
tra to estimate how much of the flux in a given bandpass 
comes from lines. However, in many cases lines are satu- 
rated in these spectra, so we instead use unsaturated, but 
uncalibrated spectra. To do this we construct a flux distri- 
bution by interpolating between the outburst JKT points 



(Hynes et al. 1998) and the mean post-outburst photom- 
etry (Clark et al. 2000). We can then multiply such a flux 
distribution by a continuum normalised spectrum and per- 
form synthetic photometry, iterating until we have a flux 
distribution such that the synthetic photometry matches 
the observed values. This is very crude, but for estimat- 
ing the contribution of lines within a bandpass, the effect 
of the assumed flux distribution is only to introduce a 



wavelength dependent weighting within the bandpass, so 
the results are weakly sensitive to the flux distribution as- 
sumed. The contribution from the lines can be estimated 
by comparing synthetic photometry performed on the flux 
distribution alone with that performed on the flux distri- 
bution multiplied by the normalised spectrum. The dif- 
ferences in magnitudes which we infer in outburst (1998 
April 4) are AB ~ 0.6, AV ~ 0.6 and AR ~ 0.9, and post- 
outburst AB ~ 0.08, AV ~ 0.10 and AR ~ 0.35. While 
the quiescent values are fairly small, and hence the qui- 
escent photometry is dominated by continuum, the out- 
burst values can be large; the R band flux, for example, is 
about 50 percent line emission, due to the strong Ha line. 
These values are similar to those quoted by Barsukova et 
al. ( p002| ). 



7. Spectral line behaviour 

The spectrum of CI Cam has long been known to be 
rich in emission lines. In outburst, this is even more dra- 
matic, with most lines showing larger equivalent widths, 
and Balmer and He 1 lines completely dominating their 
regions of the spectrum (see Fig. p]). A similarly rich out- 



burst spectrum was seen in the infrared (Clark et al. 1999). 
We now consider the evolution of selected line strengths 
and profiles before, during and after the outburst. We use 
all the data described earlier togethe r with some measure- 
ments repor ted by Wagner et al. ( 1998 ) and Orlandini 
et al. ( 2000 ). Equivalent width lightcurves for a repre- 
sentative range of lines are collated in Fig. ||; these are 
discussed in turn in the following subsections. No statis- 
tical errors have been indicated as the dominant uncer- 
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tainty is systematic due to contamination from other lines 
and the uncertain continuum level in such a rich spec- 
trum. Consequently the best indicator of the reliability of 
a lightcurve is the scatter of the plotted points. 



7.1. Hydrogen and helium lines 

During outburst Hi, Hel, and Hen lines all brighten 
dramatically. The increase is much more dramatic than 
in the continuum, so equivalent widths rise. The rising 
phase is not observed, except possibly in Hen 4686 A 
where the earliest observations suggest that the equiva- 
lent width peaked around 5 days after the X-ray outburst. 
Comparison with photometry obtained at this time indi- 
cates that this arises because during this period the He II 
flux was approximately constant while the continuum de- 
cayed. For the Hi and Hel lines, the fact that these lines 
are already decaying within two days of the X-ray out- 
burst indicates that they originate in a region that can 
respond to the outburst in less than two days. 

There are clearly major differences between the Hi 
and He I behaviours. All the Balmer lines appear to re- 
turn to approximately their pre-outburst levels within 
~ 30 days, but the Hel lines appear to drop to a factor 
of 2-10 below this. This effect has been noted previously 
by Barsukova et al. (1998), Orlandini et al. ( 200C| ) and 
Jaschek & Andrillat (2000) and is clearly real as can be 
seen in Fig. |l] where spectra before and after the outburst 
are compared. Further, there appear to be systematic dif- 
ferences between He I lines which are only manifested after 
the outburst: during outburst the Hel EW ratios are sim- 
ilar to those before the outburst, but afterwards they can 
differ, with Hel 7065 A closest to the pre-outburst level 
and He I 4713 A furthest below it. These differences are 
large, corresponding to a change by a factor of three in 
the 4713 : 6678 ratio for example, so are not simply due 
to changes in the continuum shape. This change in line 
ratios suggests changes in the physical conditions in the 
emitting gas; temperatures and/or densities; rather than 
abundance changes (e.g. due to ejection of material as 
suggested by Orlandini et al. |2000| ) . An increase in the 
6678 : 4471 ratio would be expected fro m a d ecrease in 
temperature, for example (e.g. Osterbrock 1989 ). Changes 
in the optical depth could also change ratios, and an in- 
crease in the 7065 : 4471 ratio could result from an increase 
in the optical depth (e.g. Osterbrock 1989] ). It is unfortu- 
nately impossible to be more quantitative for such opti- 
cally thick NLTE lines. Whatever changes are involved, 
they should leave the Balmer lines essentially unaffected. 

We show a selection of Hi, Hel, and Hen line pro- 
files from outburst WHT/UES data in Fig. ||. Robinson 
et al. ( |2002 ) have already discussed the profiles seen in 
outburst. The WHT/UES data were obtain ed a few days 
earlier than those of Robinson et al. ( 2002 ) but are quite 
similar. They argue that the line profiles, or at least the 
wings, are kinematic in origin and hence that emission 
seen to ~ 2500 km s -1 indicates the velocity of the hy- 
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Fig. 6. Profiles of the stronger hydrogen and helium lines 
in outburst based on WHT/UES data from 1998 April 9. 
The profiles have been rebinned by a factor of four for 
clarity. Regions contaminated by other lines have been 
omitted. 



drogen and helium emitting material outflowing from the 
sgB[e] star. There are alternative possibilities to be con- 
sidered, however. The broad wings to the lines are most 
prominent during outburst so may be associated with ma- 
terial ejected from the accreted compact object rather 
than from the sgB[e] star. Given that these lines are likely 
formed in regions of high optical depth (see Sect. 7.4), in- 
coherent Thomson scattering may also broaden the lines 
significantly; t his ha s been proposed for other sgB[e] stars 
(Zickgraf et al. 1986 ). Scattering would be expected to be 
symmetric, but the observed profiles can be interpreted 
as the sum of a broad blue shifted compone nt an d a nar- 
rower component at rest (Robinson et al. 2002 ). These 
authors obtained a satisfactory fit to the profiles with 
a double Gaussian, with the a narrow rest component 
(FWHM 50-85 km s _1 ) and a broad component (FWHM 
~ 160 km s -1 ) blue-shifted by ~ 60kms _1 . We find sim- 
ilar results from our WHT/UES spectra, although we do 
not obtain a satisfactory fit with a double Gaussian model 
to the He I lines. This is likely due to the inadequacy 
of a Gaussian fit rather than an intrinsic asymmetry. A 
double Voigt profile fit, for example, gives extra freedom 
while still using symmetric components and works much 
better. There is some sensitivity of the fit parameters to 
whether a Gaussian or Voigt profile is used, even when 
both fit well. For all fits to all lines we obtain blue shifts 
of 40-llOkms- 1 , and FWHM of 60-90 kms" 1 and 140- 
230 km s -1 for the narrow and broad components respec- 
tively, consistent with the estimates of Robinson et al. 
( |2002D . Only the blue shift, 40-llOkms- 1 , needs to be 
kinematic and we then have much lower velocities in the 
hydrogen and helium emission regions, more comparable 
to the velocities inferred for metallic lines. It remains pos- 
sible that the width of the broad components is kinematic 
however. 
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Fig. 5. Equivalent width evolution of various lines during and after the outburst. All equivalent widths have been 
normalised to pre-outburst levels, shown as a dashed line, except where noted. All points correspond to individual 
nights except for the 1998 Sept-Dec FLWO data for which monthly averages have been plotted as the evolution is 
slow by this time. The zero-point of time corresponds to the peak of the X-ray outburst, a) Balmer line evolution. Ha 
is not plotted because it is so strong as to be unreliable in many spectra, b) Hel evolution. Lines at 4921 A and 5015 A 
are omitted as they are blended with strong Fen lines in most spectra, c) The Hen 4686 A evolution, d) Evolution of 
5650-5700 A blend. This is likely dominated by Nil 5676,5686 at the peak, e) Evolution of Fen lines. The spectral lines 
used were: AA4233, 4303, 4351, 4385, 4416 (multiplet 27); AA6432, 6516 (multiplet 40); AA5197, 5234, 5275, 5316, 5325, 
5425 (multiplet 49): and AA6238, 6248, 6416, 6456 (multiplet 74). For each multiplet the sum of equivalent widths 
has been taken. For multiplets 40, 49 and 74 this was normalised to the equivalent sum before the outburst. These 
multiplets are then in approximate agreement by the end of the period observed (around day 1000), so multiplet 27, 
for which no reliable pre-outburst measurements are possible, was also arbitrarily normalised to agree with them. The 
rise between the first few observations is definitely real, being seen in all 12 lines for which measurements could be 
made, f) Evolution of [Nil] 5755 A (filled circles). 
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During the outburst, large changes are seen in these 
line profiles. This is illustrated for Ha and Hei 6678 A 
in Fig. 0. At the earliest epoch (1998 April 3) both lines 
are extremely broad and show structure not discernible 
at later times. Emission is clearly seen up to 2500 km s" 1 
in both lines, and possibly extends to ~ 5000 km s^ 1 in 
the blue wing of Ha, although this is contaminated by 
other lines. The Ha profile in fact looks remarkably like 
that of Hen S134, a near pole-on sgB[e] star in the LMC 
(Zickgraf et al. 1986). Both show a low- velocity blueshiftcd 
'notch', although the overall width of the Ha line is much 
larger for CI Cam. For Hen S134, Zickgraf et al. ( |l986|) 
suggested that the notch was due to an unresolved ab- 
sorption feature. That might also be true for CI Cam, but 
the He I profile suggests an alternative explanation. The 
latter shows inflections in both the blue and red wings. In 
fact this line clearly suggests two distinct components: a 
narrow line at rest and a broad blue-shifted line; the Ha 
profile may arise the same way, but with a broader rest 
component. This is of course exactly the decomposition 
that was applied successfully to later spectra and suggests 
that rather than using such a two-component fit as a con- 
venient parameterisation of an asymmetric profile, we can 
actually take it more literally: there really are two distinct 
emission regions. Furthermore, fitting the same two com- 
ponent model to the 1998 April 3 Hei profile gives very 
similar results to the later spectra with a component sepa- 
ration of ~ 50kms _1 and a narrow component FWHM of 
~ 75kms _1 . The main difference is that the broad com- 
ponent in the earlier spectrum is much broader with a 
FWHM ~ 750kms~ 1 . In fact, there is a hint that the 
narrow component of Ha on April 3 may itself be more 
complex. The distortion of the peak may indicate that it 
contains two unresolved components. This structure is re- 
peatable over all of the individual spectra, obtained on 
this night, and the line peak was not saturated. 

In quiescence, the weaker hydrogen and helium lines 
become narrower, approaching the width of the iron lines 
(see Figs. These lines can never be expected to be 

as narrow as the iron lines as the thermal widths will be 



larger for lighter elements. Robinson et al. ( 2002 ) observed 
broadening of 3.1kms _1 for the Fen lines, and also esti- 
mated a temperature of the iron emission region of 8000 K, 
corresponding to a thermal component of broadening of 
1.9 kms -1 . The latter will correspond to ~ 14kms _1 
for hydrogen lines. Fig. ^ shows the observed H(5 pro- 
file (which should be less affected by optical depth effects 
than stronger lines) from 2000 December 1, well after the 
Hi lines had stabilised after the outburst. We have con- 
structed a simple model profile by taking a square topped 
profile extending to ±32kms~ 1 (the model Robinson et 
al. |2002 used for the metallic profiles) and broadened it 
by ~ 14kms _1 . This model profile is somewhat narrower 
than the observed H5 profile, but the difference is not dra- 
matic, so it is plausible that the underlying Hi kinematics 
are similar to the Fe II lines. Some extra broadening may 
arise from a higher temperature in the hydrogen emission 
region and/or radiative transfer effects. 
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Fig. 7. Changes in the profiles of Ha and He I 6678 A from 
outburst to quiescence. There is clearly a dramatic change 
in the width of the lines, and the early profiles also show 
a two component structure. 



It appears that He n 4686 A is still d etecta ble after the 
outburst, as noted by Barsukova et al. ( 2002 ). A close ex- 
amination of the pre-outburst OHP spectra reveals that 
it was also present then at a higher level, albeit still weak, 
indicating that it is not simply a remnant of the outburst. 
This is somewhat surprising and another contaminating 
line cannot be ruled out. Such a line would, however, have 
to share the property of the Hei and Nil lines of being 
weaker after the outburst than before, which is not con- 
sistent with Fen line behaviour, for example. At least one 
other sgB[e] star shows Hen emission (the luminous B0 
star Hen S134; Zickgraf et al. 1986), however, so this is 
not unprecedented, and is further evidence that CI Cam 
is among the hottest sgB[e] stars. 
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Fig. 8. Change in the width of the Hei 5015 A line be- 
tween outburst and post-outburst phases. Note the simi- 
lar behaviour of the Nil 5001 A line. The upper spectrum 
has been offset by one unit vertically. 



7.2. The Bowen blend 

A common feature of the spectra of X-ray binaries is the 
Bowen blend, a mixture predominantly of N in and C in 
lines spanning ~ 4635-4650 (McClintock, Canizares & 
Tarter 1975Q , This feature does also appear to be seen in 
CI Cam, but can only be discerned in the WHT/UES out- 
burst spectra. These data are shown in Fig. |l0[ Features 
are seen corresponding to first two of the Nm 4634, 4641, 
4642 A lines; the third may be present but unresolved from 
4641 A. These features appear broader than Fe n and Cr II 
lines in the same region and are clearly declining much 
faster. It therefore seems likely that these are the Nm 
Bowen fluorescence lines. The adjacent echelle order also 
appears to show the C in 4647, 4650, 4651 A lines which 
show similar properties. Neither set of lines appear to 
be present after the outburst, although the post-outburst 
data are of lower resolution and quality. 



7.3. The 5650-5700 A blend 

There are several lines between 5650 and 5700 A which are 
dramatically enhanced during the outburst. In many spec- 
tra these are blended with other lines in the region which 
are less variable. The blend appears to be composed of 
Nil, Sen, and Fen lines (Fig. |ll|); the variable ones are 
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Fig. 9. H(5 profile from 2000 December 1. The underly- 
ing absorption (see Sect. |^) is much broader than this, 
so provides an effectively flat background which has been 
subtracted. The histogram is the data, the smooth line 
is a square profile extending to ±32kms~ 1 broadened by 
14kms _1 . The model profile has also had instrumental 
broadening applied, although this is a smaller effect. 
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Fig. 10. Nm Bowen fluorescence lines seen in outburst. 
Three normalised spectra are plotted corresponding to the 
three nights observed. The only major changes are in the 
Bowen features which become systematically fainter with 
time. The three N in components are indicated, although 
the 4642 A line, indicated by the dashed line, is not clearly 
visible and may be weaker and/or unresolved. 



likely allowed Nil lines from multiplet 73 at 5676, 5686 A. 
The strongest predicted line from this multiplet, 5679 A, 
is not clearly seen, but may be present at a lower level and 
poorly resolved from the wing of the 5676 A line. The evo- 
lution of the total equivalent width of the blend is plotted 
in Fig. |5|i and shows a rather similar behaviour to the He 1 
lines, with a lower flux after outburst than before. 
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7.5. Metallic lines 
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Fig. 11. The 5650-5700 A blend. This appears to be a 
mix of Nil, Sen and Fen lines. In the earliest outburst 
spectra the blend is unresolved, but it is clear that the 
Nil 5676 A line dominates, more so than in the spectra 
plotted here. The Nil 5679 A line should also be present; 
it may be detected in the red wing of the 5676 A line. The 
upper spectrum has been offset by 0.25 units vertically. 



7.4. Near-IR oxygen lines 

There are a number of important emission lines of O I de- 
tected in the near-IR spectrum. These were not observed 
often enough to construct lightcurves, but they still reveal 
useful information. The strongest is Oi 8446 A. This has 
an equivalent width of 41 A before the outburst and 57 A 
after it. It is very strong but as discussed by Grandi ( 1980 ) 
there are several possible reasons for this. Our observa- 
tions yield rath er similar results to those of Chakrabarty 
& Roche ( 1997 ) for GX 1+4, so we come to similar con- 
clusions: that Oi e missio n is produced by Hi Ly/3 fluores- 
cence; Clark et al. ( |1999| ) came to the same conclusion in- 
dependently from the 1.13 and 1.32 /im infrared lines. We 
find that the Oi 7774 A line is much weaker than 8446 A 
(EW ~ 2.6 A before the outburst and 4.2 A afterwards), 
corresponding to dereddened flux ratios (7774:8446) after 
the outburst of order 0.07-0.10 accounting for uncertainty 
in reddening. This rules out production of 8446 A by re- 
combination or collisional excitation, instead favouring a 
fluorescence process. Continuum fluorescence would also 
predict Oi 7254 A and 7990 A, neither of which are de- 
tectable before or after the outburst, arguing against this 
interpretation. The most likely explanation is then fluo- 
rescence by H I Ly/3. This process requires a gas which is 
optically thick to Ha, and the ratio of 8446 A to Ha can 
be used to estimate the Ha escape probability (see Grandi 



CI Cam was dubbed 'the iron star' by Downes (|1984j) and 
this is an appropriate designation. The spectrum before, 
during, and after the outburst is rich in Fen lines, in many 
cases blended in all but the highest resolution spectra. It 
is impractical to show the evolution of all lines because of 
both the large number of lines and also blending problems. 
In Fig. |^e we compile results from four multiplets covering 
a range of wavelengths and excitations. Lines from each 
multiplet that cannot be reliably deblended or that are 
sampled by few spectra are excluded. It is clear that the 
decay from outburst is extended; it may still be continu- 
ing even in the most recent observations. This is consistent 
with the observation that the iron lines appear stronger in 
post-outburst spectra than in pre-outburst ones (if they 
are still enhanced above their 'true' quiescent level). Minor 
differences appear between multiplets, e.g. compare mul- 
tiplets 40 and 49. There is no systematic trend either with 
wavelength or with excitation, however, and the origin of 
the difference is unclear. The rise in the early stages of 
the outburst is real, as it involves a difference of a factor 
of two and is repeated across 12 lines for which measure- 
ments could be made over 1998 April 3-6. It suggests that 
the iron EW peaked around day 3-10 of the outburst. 

Iron is by no means the only metallic species in the 
spectrum. For example, Robinson et al. ( 2002j ) also iden- 
tify permitted lines of Nai, Mgi, Sin, Can, Sen, Tin, and 
Cm as well as C, N, and O lines. The behaviour of these 
lines can differ dramatically. For example Fig. |l2| shows a 
selection of Fe n and Ti II lines during and after the out- 
burst. The Till lines are clearly much more variable than 
the Fen lines and are almost undetectable in quiescence. 
This is puzzling as the Fe n lines shown are of higher exci- 
tation than the Ti II lines (5.6 eV for Fe II; 4.0 eV for Ti II.), 
and the ionisation potential of Fen, 7.9 eV, is also higher 
than that of Ti II, 6.8 eV. Possibly the Ti II lines are formed 
in a different region to Fen. Cm lines may also virtually 
disappear in quiescence, although all are quite weak even 
in outburst. 

The profiles of the m etallic lines have been charac- 
terised by Robinson et al. ( 2002 ) as a square topped profile 
extending to ±32kms~ 1 and subject to a small thermal 
broadening. They attributed this to a uniformly expand- 
ing spherical shell (Fig. |l3|a). We will discuss these lines 



198C). We estimate a dereddened flux ratio (8446:Ha) af- 
ter the outburst of 0.07-0.12 and hence an Ha escape 
probability of (1.5 — 2.5) x 10~ 4 , implying a very high 
optical depth in Ha. 



in the context of sgB[e] models in Sect. 9.1, but here will 
note a complication and an alternative explanation. The 
complication is that while a square topped profile clearly 
provided a good description of the profiles at the time of 
the Robinson et al. ( 2002[ ) observations, this is not ade- 
quate at other times. Fig. |l2| shows two segments of our 
1998 April 9 WHT/UES outburst spectrum, together with 
post-outburst counterparts. The latter are at lower reso- 
lution, but are still sufficient to show deviations from the 
flat topped profile; the Fen 4297 A line clearly shows dou- 
ble peaks, or a central depression. This is essentially the 
same as the [Fen] 4287 A line, but the latter is also dou- 
ble peaked in outburst. Other lines in this spectrum show 
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Fig. 12. a) Changes in metallic lines from outburst to qui- 
escence. Tin lines are much more variable than Fen lines 
and are virtually undetectable in quiescence. The Fe n pro- 
files change, being rounded in outburst but concave in 
quiescence (like [Fen] lines). Note that the Fen line at 
4297 A in particular does not become square topped in 
quiescence. The upper spectrum has been offset by one 
unit vertically, b) Outburst spectrum of 5415-5437 A re- 
gion. The spectrum from April 14-17 is reproduced from 
Fig. 7 of Robinson et al. ( 2002 ). Note how the concave 
tops of the Cm and Fei lines in the earliest spectrum 
later become fiat tops 



similar profiles, so this is not just due to noise. The prob- 
lem is worse than this, however, as can be seen in Fig. |l2]b. 
This shows our WHT/UES spectrum of the region used 
for Fig. 7 of Robinson et al. ( 2002) ) . Two lines in partic- 
ular, Cm 5421 A and Fei 5430 A, both of which s howed 
flat topped profiles in the data of Robinson et al. (2002) 



are actually double peaked at our earlier epoch. Thus the 
deviation from a flat topped profile cannot simply be a 
post-outburst effect, and the profiles appear to evolve from 
a double peaked to flat to double peaked form through 
and after the outburst. Consequently the flat topped pro- 
files seem to be the exception rather than the norm, and 
the variations seen suggest some deviation from spherical 
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Fig. 13. Illustrative optically thin profiles for different 
outflow geometries. The viewing angle is taken to be from 
the top of the page, i.e. pole-on. a) A spherically symmet- 
ric outflow produces a square profile extending to ±Vt crm . 
b) An equatorial section of a spherical outflow also pro- 
duces a square profile (when viewed pole-on), but nar- 
rower, extending to ±Vt e rmSin#. c) If the outflow is non 
uniform, with less emission closer to the equatorial plane, 
then a central depression is created, because there is less 
low velocity emission. 



symmetry. We suggest that the underlying symmetry is 
axial rather than spherical, but that our viewing angle is 
along the axis (i.e. we see the star pole-on). If we view 
any equatorial section of a spherical outflow pole-on then 
the profile will be rectangular, with narrower equatorial 
outflows producing narrower profiles (Fig. |l3"|b). If emis- 
sion from the equatorial outflow varies with latitude then 
variations from a flat top would be seen, and a central 
depression would imply less emission near the equatorial 
plane (Fig. [l3|c 



7.6. Forbidden lines 

A number of forbidden transitions are observed in [O i] , 
[O in] , [N n] and [Fe n] . These are a rather heterogeneous 
selection of lines spanning a range of ionisation stages 
and showing ver y dive rse behaviours. As already noted by 
Robinson et al. ( 2002 ), [N n] and [O in] lines are extremely 
narrow, more so than any permitted lines (e.g. Fig. |J). 
[Ol] and [Fen] lines, however, exhibit square-edged pro- 
files similar to the Fell lines, although typically with a 
more pronounced central depression (e.g. Fig. |l2]a). In the 
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context of the equatorial outflow geometry suggested for 
the Fe II lines, a deeper central depression would imply less 
emission close to the equatorial plane. This makes sense, 
as the density is expected to be higher there and so for- 
bidden lines may well only be seen from higher altitude, 
lower density, layers. 

The line changes with time differ between species 
too. Unfortunately, many of the lines are weak and/or 
blended with other lines, principally of Fen. The only 
line well suited to quantitative analysis is [Nil] 5755 A, 
as it is strong and unblended. The equivalent width evo- 
lution of this line is shown in Fig. ^|f. The early evolu- 
tion 10 days) is roughly consistent with a constant line 
flux, with the equivalent width increasing as the contin- 
uum fades. Comparison with outburst photometry (Clark 
et al. 2000 ), after correcting for the increased contribu- 
tion of the lines to broad band flux (Section ||) , indicates 
that this outburst line flux is comparable to that seen in 
1998 January. Barsukova et al. ( 2002 ) note that the line 
flux does increase modestly around 50-250 days after out- 
burst, then decays again. The decay can clearly be seen 
in Fig. |)f, although the rise is hard to distinguish from 
that due to the decreasing continuum flux without flux 
calibrated spectra. 

The [O in] 5007 A line also exhibits a narrow profile 
but is clearly stronger in outburst than after it (Fig. ^J). 
The 4959 A line behaves in a similar way. These lines are 
difficult to measure in most spectra as they are weak and 
blended, but it is clear that the outburst flux is much 
greater than pre-outburst and that after the outburst they 
decline, but like [N n] , they do not appear to have dropped 
to the pre-outburst level. For example, for [Oiii] 5007 A, 
the pre-outburst EW is difficult to measure but - 0.25 A. 
On 1998 April 11 it had an EW of 2.0 A and by 1998 July 
20 this had dropped to 1.1 A. As late as 2001 October 23 
it was still > 0.8 A. The decline in EW from 2000 April 
to July corresponds to a decrease of about a factor of 
three in line flux. Since the [O in] lines have very long re- 
combination times, this implies that significant collisional 
de-excitation must be occurring, possibly when the radio 
ejecta reach the line formation region. Assuming an ex- 
pansion velocity of ~ 5000 km s -1 , then a reduction of a 
factor of three within ~ 100 days implies that most of the 
[Oiii] emission originates within ~ 300 AU. 

[Fe n] lines also strengthen moderately during outburst 
(e.g. Fig. |l2|a) and [O i] lines increase enormously, being 
almost undetectable before and after the outburst. 

8. Short term spectral variability 

There have been some claims of short term variability 
in CI Cam during outburst. At X-ray energies Frontera 
et al. Q1998D found with Beppo-SAX that the 0.5- 
1.0 keV lightcurve showed significant variations on ~ 100 s 
timescales on 1998 April 9-10. They found no significant 
variation from a smooth decay in the simultaneous 1.5- 
10 keV data however, nor were any variations seen on April 
3. Ueda et al. (|1998|) examined ASCA data from 1998 



April 3-4. They also found that there was no variability 
above 1 keV, but that soft flares were seen below 1 keV on 



timescales of a few hours. Belloni et al. (1999) examined 
RXTE data spanning April 1-9 and found no evidence 
for any variation other than a smooth decay at any time. 
The short term X-ray variability thus appears confined to 
flares in the soft (^ 1 keV) band. 

Optical variability appears sporadic too. Frontera et 



al. (1998) also found evidence for 0.3 mag optical flickering 
on hour timescales on April 6, but not in later observa- 



tions spanning April 10-26. Clark et al. ( 2000 ) examined 
photometry from April 13 and 19 and found no variations 
with amplitude greater than 1 percent. 

The McDonald spectra obtained during the decline 
from outburst were taken as a series of short exposures 
to facilitate studies of line variability. Because a relatively 
narrow slit was used, however, there are significant slit 
losses. Fortunately, as noted above (Clark et al. 2000 ), 
there was no detectable short-term variability in simulta- 
neous R band photometry, so we chose to normalise each 
spectrum before extracting emission line lightcurves; these 
lightcurves are effectively equivalent widths. We find no 
evidence for variability on timescales from 3 min to ~ 2 hr 
in Ha or Fell lines. Formally we find rms scatters in the 
lightcurves of 1.1-1.6 percent in Ha and 0.8-1.5 percent for 
a composite of several iron lines. The one line that may ex- 
hibit short term variability is Hei 6678 A. This shows 1.9, 
4.9 and 3.5 percent rms variations on 1998 April 18, 19 and 
20 respectively. On the latter two nights there is clearly 
a systematic variation. To be sure this is not an artifact 
we have renormalised the continuum using just the 6648- 
6663 and 6693-6708 A regions and recentered the line in 
each spectrum by cross-correlating the line profiles with 
an average. None of these changes affect the result: this 
line appears to show real variability. The case is most per- 
suasive on the last night. There is a clear overall rise in the 
lincstrength of ~ 10 percent over < 3 hours, shown in Fig. 
|l4| . This is in the opposite direction to the overall decline 
in the linestrength and much larger than expected from 
the changes in the continuum strength alone (e-folding 
decay time ~ 24 days at this time; Clark et al. 2000). We 
have also constructed rms spectra for each night and com- 
pared them with the average spectra. These suggest the 
same conclusion as the lightcurves; the only feature that 
seems to show significant variability is the He I line on the 
second and third nights. 

We also examined the FLWO data from 1998 April 
3 in the same way. The rms spectrum shows no variable 
features other than Telluric bands. The line strengths for 
Ha and Hei 6678 A show an rms scatter of 1.8 percent 
and 1.2 percent respectively with no systematic trend. We 
believe this represents a null detection of variability on 
this night. 



9. Discussion 
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Fig. 14. Change in the strength of the Hei 6678 A line in 
McDonald spectra from 1998 April 20. Spectra were nor- 
malised to a flat continuum which was then subtracted. 
Integrated line fluxes are plotted as fractional variations 
about the mean. Error bars are formal errors obtained 
in extracting one-dimensional spectra and do not include 
additional uncertainties that may be introduced by nor- 
malisation of the spectra. 

9.1. CI Cam as an sgB[e] star 

It is clear that the X-ray outburst was associated with dra- 
matic changes in the spectrum of CI Cam. Having iden- 
tified it as an sgB[e] star, however, it will prove useful to 
first compare it with the small number of other stars of 
this class and see how much of its 'unusual' behaviour, 
is typical of these objects, and what requires additional 
input from a comp act ob ject. 

Robinson et al. ( 2002 ) infer a predominantly spherical 
geometry for the outflow from CI Cam. Several compo- 
nents are suggested. A cool, low velocity (32kms _1 ) wind 
gives rise to the many iron lines and the square profiles 
of these lines suggested that this component is spherical. 
Multiple higher velocity components (> 1000 krns" 1 ) are 
suggested by hydrogen and helium lines, and by the UV 
resonance lines. Thi s pictu re has several problems as noted 
by Robinson et al. (2002), principally i) it is unclear how 
several near-spherical winds can co-exist with very differ- 
ent velocities and ionisations; and ii) the quiescent X-ray 
luminosity is very low for a compact object continually 
burrowing through a dense, spherical outflow. 

This picture is also rather at odds with existing mod- 
els of sgB[e] stars (Zickgraf et al. 1985 1986; Oudmaijer 
et al. 1998[ ). The structure inferred from other systems is 
of a two-component wind which is far from spherical. A 
rarefied, hot, high velocity wind dominates in the polar 
regions and is responsible for the UV lines often seen with 
P Cygni profiles. This wind is identical to the winds of nor- 
mal early- type supergiants. Other lines are attributed to a 
cooler, denser wind concentrated in the equatorial regions 
that is not present in normal supergiants. The nature and 
origin of this equatorial material remain unclear, and it is 



probably quite different to the discs around classical Be 
stars. Can CI Cam be fitted into this framework? 

The interpretation of the UV lines seen in CI Cam 
is most straightforward. These show broad P Cygni lines 
typical of supergiant winds. They can be associated with 
the hot, high velocity, polar outflow, and we essentially 



agree with Robinson et al. (2002) on the interpretation of 
these lines. 

The hydrogen and helium lines are more problem- 
atic, however. In outburst, these lines are broad and 
asymmetric, with blue-shifted emission extending to J> 
2500kms _1 , but with no d etectable absorption compo- 
nents. Robinson et al. ( 2002] ) associate these lines with a 
high- velocity, weakly collimated outflow from the sgB[e] 
star. It is unclear where this could be situated, however: 
if it originated from the same region responsible for the 
UV lines we might expect P Cygni line profiles, as seen in 
HD 87643 (Oudmaijer et al. |1998| ). The latter object may 
be rather different to other sgB [e] stars with a higher H I 
column density in the polar wind. In most sgB[e] stars, 
hydrogen and helium emission is instead attributed to the 
equatorial component, and hence relatively low velocities 
should be involved, as is seen in the iron lines. In fact, 
a clue that the hydrogen, helium and iron lines are all 
associated with the same region, or regions, is provided 
by similarities of the line profiles. Robinson et al. (2002) 
associate Fell lines with a very different region to the 
hydrogen and helium lines, but they also note that 'The 
stronger [iron] lines have more rounded profiles, becoming 
almost Gaussian in shape as the lines become stronger, 
and the very strongest lines have an extended blue wing. ' 
Another way to put this is that the strongest lines be- 
come more akin to hydrogen and helium lines. While the 
rounded profiles could indicate radiative transfer effects, 
the presence of asymmetry in the Fe n lines suggests that 
the same rest and blue-shifted components are present in 
the iron lines as are seen in hydrogen and helium pro- 
files. If the strongest, optically thickest iron lines look like 
hydrogen lines, then the corollary is that the weakest, op- 
tically thinnest hydrogen and helium lines should look like 
iron line s. T his does seem to be the case, as demonstrated 
in Sect. [Ft] and Fig. | 

The decomposition of the hydrogen and helium lines 
into two clear components suggests that two distinct re- 
gions are involved during outburst. The narrow rest com- 
ponent dominates in the iron lines and is also important 
in the hydrogen and helium lines. A second broad, blue- 
shifted component is present in hydrogen and helium lines, 
and weakly detectable in the strongest iron lines. The nar- 
rower, rest component is likely the same component as 
seen in quiescence, and in other sgB[e] stars, originating 
from the equatorial outflow. The profiles seen in Hi and 
He I are very different to the Fen lines, but the former 
are broadened by larger thermal widths and inco herent 
Thomson scattering, as suggested by Zickgraf et al. ( 1986 ) 
for other sgB[e] stars. The strong, broad, blue-shifted com- 
ponent is a feature of the outburst and may not be associ- 
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ated with the outflow from the sgB[e] star at all. We will 
discuss this possibility further in the following section. 

Associating the iron lines with the equatorial material 
conflicts with the interpretation of their profiles offered by 
Robinson et al. (2002). They argue that the rectangular 
profiles allow little deviation from spherical symmetry and 
no significant rotational velocity. As we have discussed in 
Sect. 



7.5 



however, large deviations from spherical geom- 
etry are possible in one case: that an equatorial outflow 
is viewed pole-on. A pole-on equatorial outflow can also 
explain the apparent lack of rotational velocities, as the 
projected rotational velocity could be very low, and the 
intrinsic rotational velocity in regions of the disc where 
Fe II can arise is in any case likely to be very low. 

As a further argument for the presen ce of a spherical 
low velocity wind, Robinson et al. ( |2002 ) claim on the ba- 
sis of the scaling relation of Bjorkman ( 1998 ) that dust 
can condense in this wind and that an equatorial outflow 
is therefore not required for the production of dust. The 
scaling relation of Bjorkman ( 1998 ) is based on the ar- 
gument that the density of wind compressed outflows in 
sgB[e] stars at radii where the wind is cool enough for 
dust to condense is ~equal to those of post-AGB stars 
(which do form dust) and therefore that it is possible for 
dust to condense in such outflows. However, the scaling 
relation does not allow for the fact that the winds in post- 
AGB stars are significantly enhanced in species such as 
Si and C (via dredge up) which form dust up to a fac- 
tor of ^10 times more easily relative to the winds from 
sgB[e] stars. Additionally, the mass loss for the sgB[e] star 
is calculated from the base density of the outflow at the 
surface of the star and the terminal velocity of the outflow, 
rather than the velocity of the outflow at the surface of 
the star. This will also serve to overestimate the density of 
the outflow at the dust condensation radius. Combining 
the 2 arguments we estimate that the density is likely to 
be 2-3 orders of magnitude too small to allow dust to 
condense (J. Porter, priv. comm.). While this does not a 
priori argue against the presence of a quasi-spherical wind 
resulting in the Fe n line profiles it does demonstrate that 
dust cannot condense in it and that a region of higher 
density - whether an equatorial outflow or dense ejecta 
- is required for condensation. Robinson et al. ( 2002 ) ad- 
ditionally argue that the presence of a resolved, circular 
dust shell (Traub, Millan-Gabet & Garcia 199S| ) is consis- 
tent with their interpretation, but we note that if a disc is 
viewed pole-on then dust in the outer regions of the disc 
will produce the same apparent structure. 

Forbidden lines are a rather heterogeneous mixture, 
but broadly divide into two types. Narrow lines ([Nil], 
[O in] ) are formed in a region of very low bulk, turbulent 
and/or thermal velocities. This probably places them out- 
side the two component outflow from the sgB[e] star, and 
they may be a remnant of an earlier p hase of evolution 
of one of the stars (c.f. Oudmaijer et al. 1998 ). The [Fen] 
lines show profiles similar to the Fe n lines, but with a more 
pronounced central depression. As discussed in Sect. |7l| 
this suggests that these lines may be formed in the upper 



layers of the same equatorial material responsible for Fen; 
[O i] lines are likely formed in the same region. 

To summarise, by analogy to other known sgB[e] stars 
we suggest that when the system is not in outburst, 
the strongest optical permitted lines (Hi, Hei, Fen, and 
others) and some forbidden lines ( [Fe n] , [O i] ) arise in 
an equatorially concentrated region, whether a quasi- 
Keplerian disc (e.g. Okazaki 2001), a hybrid disc with 
a slow radial expansion velocity, or even in a radiation 
driven disc wind such as is suggested in the Galactic sgB[e] 
star HD 87643 (Oudmaijer et al. |1998D . UV P Cygni lines 
originate from a hotter, higher velocity polar outflow of 
much lower density. The whole structure is viewed close 
to pole-on, hence disc lines have rather low velocities as 
both rotational and expansion velocities would be mainly 
directed in the equatorial plane. Finally the narrow forbid- 
den lines ([Nil], [Oin]) come from a much more extended 
and near stationary region. 

Our suggested model does depend on viewing the 
source near pole-on. There are several independent lines of 
evidence to support this, suggesting that CI Cam is seen 
at a lower inclination than many sgB[e] stars. 

1. Most sgB[e] stars show low velocity blueshifted ab- 
sorption in their Balmer line profiles (Zickgraf et al. 



1985,1986; Oudmaijer et al. 1998). This can be inter 



preted as due to absorption by the slow equatorial wind 
component. No such absorption is seen in CI Cam sug- 
gesting that the inclination is low enough that our line 
of sight does not pass through the slow wind. 
sgB[e] stars are expected to show significant po- 
larisation changes through Ha. The observations of 
HD 87643 are qualitatively similar to what is ex pecte d 
for a rotating, expanding disc (Oudmaijer et al. 1998 ). 
In CI Cam no polarisation changes are seen through 
the lines (Ikeda et al. 2000 ). In the models of Wood et 
al. ( 1993| ) the polarisation decreases with inclination, 
so a lack of line polarisation also supports a low inclina- 
tion. In addition, the outburst observations of Ikeda et 



al. (2000) showed a continuum polarisation consistent 
with surrounding field stars, and hence with purely in- 
terstellar polarisation. This is also as expected for a 
nearly pole-on disc. 
3. We have argued above that dust cannot form in a 
spherical shell and that it must be associated with the 
outer disc. The near perfectly circular resolved infrared 
image then suggests that this disc must be viewed pole- 



It therefore seems plausible that CI Cam is a pole-on 
sgB[e] star. Such a model is consistent with the properties 
of both the optical and UV spectra. 

9.2. The nature of the outburst 

We now move on to discuss the role of the compact ob- 
ject. Its nature obviously has some bearing on this discus- 
sion. The X-ray outburst was very short, much shorter 
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than typical of soft X-ray transients (SXTs). This, to- 
gether with the radio emission from an expanding rem- 
nant (Mioduszewski et al. 2002, in preparation) suggests 
some kind of explosive event. Orlandini et al. (2000) sug- 
gest a thermonuclear runaway on the surface of a white 
dwarf is responsible, but this seems hard to reconcile with 



the high inferred X-ray luminosity (Robinson et al. 2002) 
and the 7-ray detection by CGRO/B ATSE (Belloni et al. 



199E). It seems more likely that the outburst involved a 



brief burst of supercritical accretion onto a neutron star 
or black hole, resulting in ejection of much of the accreted 
material. This accounts for the observed radio ejecta and 
possibly the broad components of the hydrogen and he- 
lium lines. A supercritical accretion model could also in- 
volve a large enough optical depth of scattering material 
around the X-ray source to smear out any short timescale 
variability. 

One possible cause for such a large burst of accre- 
tion would be the passage of the compact object through 
the equatorial plane. We can attempt to estimate the ac- 
cretion rate, although this will be very uncertain as the 
orbital parameters are unknown and the physical con- 
ditions in sgB[e] star outflows are not well understood. 
We assume stellar parameters typical of hot, luminous 
sgB[e] stars (Zickgraf et al. 1986) and adopt a representa- 
tive 30 year orbit of a 10 M Q black hole around a 60 M Q 
sgB[e] star with an eccentricity of e ~ 0.9. The perias- 
tron distance is then about 10-20 R s tar and the perias- 
tron velocity 170 km s -1 . In the model of Oudmaijer et 



al. ( J1998[ ) for HD 87643, the equatorial density at 10 R star 
is ~ 2.5 x 10~ 12 gcm~ 3 . This model may not be correct, 
and CI Cam is probably more massive and luminous than 
HD 87643, but this at least indicates the kind of equa- 
torial density which is considered plausible. The velocity 
of the compact object relative to the equatorial material 
will be dominated by its orbital motion, so the Bondi ac- 
cretion rate is then predicted to be ~ 6 x 10 21 gs _1 or 



~ 400MEdd (Bondi 1952 ). Given the uncertainties men- 
tioned above, this is a very approximate estimate, and 
could be at least an order of magnitude off. Even allowing 
for this large uncertainty, the accretion rate in this sce- 
nario can therefore be expected to be extremely high, due 
to the high density and low velocity of the equatorial mate- 
rial, and hence a supercritical accretion episode is possible. 
In contrast, when the compact object is out of the plane 
the density is lower by a factor of 10 4 (in the model of 
Oudmaijer et al. 1998) and the relative velocity is higher, 
since the high latitude outflow moves much faster. The 
Bondi accretion rate is then expected to be much lower, 
~ 8 x 10 -5 MEdd ■ If the accretion efficiency remained high, 
77 ~ 0.1, then CI Cam should still be a relatively bright X- 
ray source, with Lx ~ 10 35 ergs -1 . However, at these low 
accretion rates then the flow could be expected to become 
advective as proposed for other quiesc ent bl ack hole candi- 
dates (Narayan, Garcia & McClintock 2001 and references 
therein), and the accretion efficiency would then be lower, 
r\ ~ 10~ 2 -10~ 4 , imply ing L x ~ 10 32 ~10 34 ergs" 1 as ob- 
served (Robinson et al. 2002 ). This model therefore avoids 



the problem that Robinson et al. ( |2002 ) had, that a com- 
pact object burrowing through a dense spherical outflow 
should be persistently bright. 

The passage of a compact object through the equato- 
rial material, and the ejection of a lot of material from 
near it could affect the equatorial region significantly, so 
changes in the line spectrum from this region would be 
expected. We cannot offer an exact mechanism for this, 
however, and several factors may be involved: the tidal 
effect of the compact object passage; X-ray heating; and 
the interaction of the expanding radio remnant with the 
equatorial flow. The rapid response of the lines to the out- 
burst, peaking within a few days of the X-rays, constrains 
the tidal effect, as a tidally triggered disc outburst would 
be expected to proceed on a viscous timescale, if the disc is 
Kcplcrian (i.e. thin and rotationally supported, with small 
inflow or outflow velocity). This is likely to be very long, 
hundreds to thousands of days. It is far from clear that 
sgB[e] discs are Keplerian, but other timescales, e.g. the 
outflow time from the stellar surface to ~ 50 stellar radii, 
are likely to be similar. The extended decay of the iron 
lines may indeed be due to a recovery of the disc to its 
pre-outburst state on this timescale, but the rise of the 
outburst cannot be, so a tidally triggered disc outburst is 
unlikely. 

The X-ray heating effect will obviously be much faster, 
and should rise and fade on a timescale comparable to the 
X-ray outburst itself, although could be prolonged by sig- 
nificant cooling or recombination times. This is not seen, 
although a signature of X-ray heating of the disc may be 
present in the form of strong 6-7 keV line emission seen 
by SAX (Frontera et al. |l998|), ASCA (Ueda et al. ^998|) 
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and RXTE (Revnivtsev et al. [199J ; Belloni et al. 
If this line is attributed to fluorescent iron emission then 
a large, cold, optically thick surface is needed covering 
about half the sky as seen from the X-ray source (Ueda 
et al. 1998 ). This is what would be expected if the X-ray 
source were above the plane of the equatorial material. 
The dramatic enhancement of high excitation lines, such 
as He 11 and N ill during outburst also suggests that X-ray 
or extreme- UV irradiation may place a role. 

The expansion of the radio ejecta can more naturally 
produce the timescales observed than either a tidal inter- 
action or X-ray irradiation can. Assuming an expansion 
at ~ 5000 km s -1 (from Mioduszewski et al. in prepara- 
tion, assuming d ~ 5kpc), this would reach a disc radius 
of ~ 50 stellar radii (~ 2000 Rq) in about 3 days, so this is 
consistent with the relatively short outburst and the rise 
time of the Fe 11 EWs. The expanding ejecta may also ex- 
plain the broad blue-shifted component seen prominently 
in hydrogen and helium profiles during outburst, as the 
inferred expansion velocities are comparable to the max- 
imum detectable velocities in the emission lines. This is 
an appealing interpretation, as the origin of these ejecta 
is most likely the compact object, and hence an overall 
blue shift can be explained as the orbital velocity of the 
compact object at the time of ejection. The overall blue- 
shift seen, 50-100 km s -1 , is reasonable compared to our 
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estimate of the periastron velocity ~ 170 km s -1 . It is also 
sensible that we see a blue-shift, rather than red, as X- 
ray spectra taken in the outburst decay show a strong 
fluorescent iron line, consistent with reflection from the 
equatorial material, and do not show strong absorption; 
both indicate that the compact object was likely between 
us and the equatorial plane at the time, and hence should 
have been moving toward us. 

9.3. The origin of the optical continuum 

The dramatic increase in the optical luminosity of CI Cam 
during the outburst is also intriguing. The optical rise was 
not observed, but if it was later than the X-rays then the 
rise time was fast, Sa 2 days. While large optical brighten- 
ing is typical of SXTs, it is actually rather surprising in a 
system containing a hot, luminous sgB[e] star. The bright- 
ening is such that whatever provides the additional source 
of light during outburst must significantly outshine the 
sgB[e] star itsel f. Rel ative to the mean pre-outburst level 
(Bcrgner et al. 1995 ) the brightest magnitudes re porte d 
are brighter by 1.9 mag (U, Apr 3.9; H ynes et al. |l998| ), 
2.1 mag (B, Apr 3 .1; Ga rcia et al. |1998| ), 2.3mag (V, Apr 

98), 3.5 mag (R, Apr 2 .1; R obinson 

19981 ), 2.4mag(7, Apr 3.8; Clark et al. |2000| ). These 
observations are not simultaneous and are not intended 
to indicate the spectrum; merely to indicate that the out- 
burst was dramatic throughout the optical region. Since 
the peak of the outburst was missed in the optical these 
are actually lower limits on the outburst amplitude. While 
some of the brightening is due to stronger line emission, 
not all can be; removing the line emission from our early 
spectra only increases B and V by ~ 0.6 mag, so much 
of the brightening has to be enhanced continuum. Fig. ^ 
indicates that on 1998 April 3, the excess continuum was 
2-5 x brighter than the sgB[e] star, and redder than it. 
Where can this extra continuum come from? Some possi- 
bilities which we discuss below are: 

1. Heating of the sgB[e] star during the outburst. 

2. Heating of the dusty envelope around the system. 

3. Emission from an accretion disc around the compact 
object. 

4. Enhanced emission from the circumstellar material 
around the sgB[e] star. 

5. Emission related to the observed radio ejection. 

To help constrain the possibilities it is useful to ex- 
amine the spectrum of the excess light. Fig. [l5| shows the 
difference between the early outburst and post-outburst 
spectra from Fig. || Close examination of the outburst 
spectrum, and comparison with high resolution spectra 
taken later in outburst, indicates that although much of 
the spectrum is dominated by emission lines, the lowest 
points do appear to indicate the continuum. The differ- 
ence spectrum is obviously rather uncertain, but many 
possibilities can be ruled out. If it is a power-law, then 
it is relatively flat in /„. The power-law model plotted is 



/„ oc v 1 ! 2 . A completely flat spectrum would also be pos- 
sible, but it could not be much redder or bluer than these 
possibilities. A hot black body is ruled out, as indicated 
by the fact that the difference spectrum is redder than 
the sgB[e] spectrum. A 10 000 K black body is shown; a 
temperature much cooler than this would also not be con- 
sistent with the data. Similar conclusions can be drawn 
by examining the change in B — V during the outburst, 
from ~ 0.8 mag to Ss 1.2 mag. It is obviously difficult to 
be precise given the uncertain reddening and large line 
contribution, but a very hot or very cold black body, or a 
very red or very blue power-law can be ruled out. 

This colour information can immediately be used to 
constrain the origin of the excess light. Interpretation 1 
can be ruled out as heating of the sgB[e] star could not 
make the spectrum redder even if a 5 x increase in bright- 
ness were possible. In this case, the excess should be ex- 
tremely blue. Equally, interpretation 2 can be ruled out 
as the excess light is not cool enough. Heated dust would 
dissociate for T ^ 2000 K, but the observed spectrum and 
change in B — V are clearly not consistent with such a low 
temperature. 

The brightness of the outburst argues against interpre- 
tation 3. After correcting for line emission (AV^ ~ 0.6 mag) 
and extinction (Ay ~ 4 mag), the peak dereddened con- 
tinuum magnitude is V 5.8 and the average quiescent 
value is (V) ~ 7.7 mag; hence the additional continuum 
source has V ^ 6.0. At a distance of ~ 5kpc this corre- 
sponds to an absolute magnitude of My —7.5. This is 
much brighter than any low mass X-ray binary (LMXB) 
(—5 < My < +5; van Paradijs 1994), i.e. much brighter 
than other observed accretion discs around black holes 
or neutron stars. If the additional light were an accre- 
tion disc it would then have to be more luminous than 
those in LMXBs. Since the spectrum indicates that it is 
not extremely hot, it would have to be large. As the size 
of LMXB accretion discs is limited by tidal truncation, 
a much larger disc in CI Cam is certainly possible. If the 
density in the disc were comparable to the disc densities in 
LMXBs, however, such a large disc would provide enough 
mass to sustain a much longer X-ray outburst than seen, 
and a decay timescale of months, more like that seen in 
SXTs, would be expected. Indeed Robinson et al. ( 2002 ) 
have argued that the very rapid decay timescale in CI Cam 
indicates that the accretion disc, if present, must be very 
small. This is a very different environment from that in an 
LMXB, however, and the disc could be stabilised against 
thermal instability by irradiation from the supergiant at 
much lower densities than required in an LMXB, resulting 
in a lower disc mass. The brightness remains a problem, 
however, as the disc luminosity would be comparable to 
that of the supergiant and greater than that of the X- 
ray source at the same time, and so reprocessing of light 
from either is unlikely to dominate. It therefore seems un- 
likely that an accretion disc around the compact object 
can dominate the optical brightening in outburst. 

The relatively low temperature and large area inferred 
for the excess light, if it originates in black body emission, 
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can more plausibly be associated with the equatorial ma- 
terial around the sgB[e] star than an accretion disc around 
the compact object. For a temperature of ~ 10 000K, 
the extra emission requires an area ~ 40 x that of the 
sgB[e] star. Both the temperature and size are plausible 
for the equatorial material around a sgB[e] star. In this 
case the timescale arguments are similar to those for the 
lines. The optical decay appears slower than the X-rays 
(e.g. Clark et al. 2000). The peak optical luminosity of 
the continuum source is also large; for T ~ 10 000 K and 
an area ~ 40 x that of the star, the luminosity was at least 
- 2x 1O 5 L - 8x 10 38 ergs _1 (assuming d - 5kpc). This 
is actually larger than the peak X-ray luminosity at this 
distance, Lx ~ 3 x 10 38 ergs _1 , and by the time of the 
optical observations, the X-ray flux had already d roppe d 
by more than an order of magnitude (Belloni et al. 1999|) . 
Both the decay timescale and optical luminosity render it 
unlikely that the continuum brightening can come from 
X-ray reprocessing. The timescales would be more plausi- 
ble for an interaction of the equatorial material with the 
radio ejecta; assuming an expansion speed ~ 5000kms _1 , 
this area will be covered in less than half a day, so a rapid 
optical brightening is possible in this way, i.e. powered by 
the kinetic energy of the ejecta. 

The remaining possibility is that there is a consider- 
able direct optical contribution from the ejecta from the 
explosion. The similarity of timescales of the optical and 
radio decay (Clark et al. 2000) argue for this interpre- 
tation. Optical synchrotron can probably be ruled out. 
The dereddened magnitude of the additional continuum 
component of V Ss 6 corresponds to a flux of <; 15 Jy 
early in the outburst. Contemporaneous radio observa- 
tions were much weaker than this, with the spectral energy 
distribution peaking at 650 mJy at 8.4 GHz and de creasin g 
at higher frequencies (Hjcllming & Mioduszewski |l998b|) . 
Optical synchrotron emission thus seems unlikely without 
a very unusual electron energy distribution. However, if 
the accretion rate was highly supercritical, we might ex- 
pect much of the energy released to be reprocessed by 
a scattering envelope to produce a bright optical source 
(Shakura & Sunyaev 1973). Such a supercritical accretion 
scenario has been invoked for a number of X-ray binaries. 
SS 433 i s like ly to a persistently supercritical source (e.g. 
Fabrika 1997 and other authors; see also Okuda & Fujita 
200C and references therein). V4641 Sgr, which also had 
an extremely short but luminous X-ray outburst has been 
suggested to have undergone a transient burst of supercrit- 
ical accretion (Revnivtsev et al. 2002). Shakura & Sunyaev 
(1973) predict that the optical spectrum in the supercriti- 
cal regime should appear as a power-law, F y cx v 1 / 2 , satu- 
rated with broad emission lines. We have already discussed 
the origin of the broad emission lines that we see and sug- 
gested that these may be related to outflowing material 
from supercritical accretion. The outburst optical contin- 
uum emission is clearly redder than the underlying star, 
and the excess continuum emission is consistent with a 
F v oc v 1 / 2 power law (Fig. [l5|). This therefore does seem 
a plausible interpretation. If this supercritical accretion 
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Fig. 15. Spectrum of the extra light seen in outburst, ob- 
tained by taking the difference between spectra obtained 
on 1998 April 3 and 1998 October 29. This has been dered- 



dened using the Fitzpatrick (1999) extinction curve as- 
suming Ay = 4.0. Two possible models for the continuum 
spectrum are also plotted; the power law corresponds to 
F v oc v 1 ' 2 . 



scenario is correct, then the observations suggest that the 
accretion rate was extremely high and hence the spheri- 
sation radius large; assuming a terminal velocity for the 
ejected material of ~ 5000kms _1 , and an optical luminos- 
ity ~ 10 38 ergs -1 , suggests an accret ion rate of M ~ 10 2 - 
10 4 M cr it (Shakura & Sunyaev 1972), consistent with our 
estimate above of the possible periastron accretion rate. 

All of these calculations are extremely simplistic. Our 
aim is not to present a detailed spectral model, but rather 
to test which explanations of the outburst emission are 
plausible. X-ray heating, whether of the sgB[e] star it- 
self, the equatorial material, an accretion disc around the 
compact object or extended dust, is not consistent with 
observations. Heating of the equatorial material by the 
interaction with the expanding radio remnant, or direct 
emission from these ejecta, remain possibilities. The lat- 
ter has the advantage that invoking a supercritical accre- 
tion regime explains both the optical continuum emission 
and the broad, blue-shifted emission components, so is our 
favoured interpretation, but it is likely that a combination 
of these mechanisms is involved. 



9.4. Is CI Cam an X-ray binary? 

Finally we emphasise that none of the preceding discus- 
sion depends on CI Cam being a binary system. We, and 
others, have argued that CI Cam is an sgB[e] star, and a 
compact object of some kind is clearly implicated in the 
X-ray outburst, but it is not necessary for them to be 
physically associated. Indeed, the properties of CI Cam 
in quiescence are sufficiently similar to other hot, lumi- 
nous sgB[e] stars that if it is a binary then the compact 
object's influence appears to only be important during an 



24 



R. I. Hynes et al.: Spectroscopic observations of CI Cam 



outburst. For example, it clearly does not truncate the 
disc as can happen in classical Be X-ray binaries. While 
it could be that CI Cam does have a compact object in 
a long period orbit, it is also possible that the outburst 
could have resulted from a chance encounter of an isolated 
black hole or neutron star with the circumstellar material. 
As a massive, young object, CI Cam is likely located in 
a region of recent star formation, so there are likely to 
be many stellar remnants in its proximity. While such a 
chance encounter still seems improbable, it cannot be dis- 
proved based on the statistics of a single event. So far we 
have no conclusi ve evi dence for the binarity of CI Cam. 
Barsukova et al. ( 2002 ) have suggested that two periodic- 
ities are present. Their short 11.7 day period is certainly 
too small to represent the orbital period in a system such 
as this as the compact object would have to be almost on 
the surface of the sgB[e] star, or even inside it. The longer 
1100 day period would be more plausible, but cannot be 
considered convincing until it is seen to repeat for multiple 
cycles. In balance, while the hypothesis that CI Cam is a 
binary seems most likely, it is not proven and the alterna- 
tive, that there was a chance encounter with an isolated 
compact object, cannot be ruled out. 

10. Summary and conclusions 

CI Cam is an sgB[e] star, and many of its characteristics, 
particularly the emission line spectrum, are representative 
of this class. Comparison with other sgB[e] stars suggests 
that CI Cam is among the hottest members of the class 
and viewed nearly pole-on. It differs from the other sgB[e] 
stars in interacting with a compact star which introduces 
an additional variable element. It seems most likely that 
the compact object is physically associated with CI Cam, 
i.e. that it is an HMXB. In this case, the compact ob- 
ject is probably in a long period, eccentric orbit. However 
we cannot rule out the possibility that this was a chance 
encounter and that the two stars are not physically asso- 
ciated. Resolution of this issue will likely involve waiting 
for a true periodicity, repeated over several cycles, or an- 
other outburst. In considering the outburst mechanism, 
there is actually little difference between a long period ec- 
centric orbit and a chance encounter, so our discussion of 
the outburst applies to both cases. 

We suggest that the majority of the optical emission 
lines originate from an equatorially concentrated outflow 
or circumstellar disc. During outburst, hydrogen and he- 
lium lines appear to have two components, a narrow rest 
component and a moderately blueshifted broad compo- 
nent. Metallic lines are mainly dominated by the narrow 
component at all times, although some asymmetry is seen 
in outburst suggesting that a broad component is present. 
The square profiles of the Fen lines can be explained by 
the equatorial outflow model, if viewed pole-on, so the nar- 
row component is likely associated with this. The broad 
component becomes weaker and narrower on the decline, 
and almost disappears in quiescence. This may come from 
material ejected in the outburst rather than from the 



equatorial outflow. Forbidden lines fall into two categories; 
[Oi] and [Fen] show similar line profiles to Fen, but with 
a lack of low velocity material. These profiles could orig- 
inate from the low density, upper layers of the equatorial 
outflow. Other forbidden lines, [Nil] and [Oni] are nar- 
rower and likely come from a much more extended region. 

The outburst mechanism remains undetermined, al- 
though the outburst was probably precipitated by the pas- 
sage of the compact object through the equatorial mate- 
rial. It is unlikely that X-ray heating of any component 
is responsible for the optical outburst. Instead the optical 
outburst is likely associated with the expanding remnant 
produced by the X-ray outburst, either through direct 
emission from the remnant or as a result of its interac- 
tion with the circumstellar material. The spectral shape 
of the outburst optical continuum, and the presence of 
broad, blue-shifted emission components, are both consis- 
tent with predictions for supercritical accretion resulting 
in ejection of much of the material (Shakura & Sunyaev 
1973), and the peak mass transfer rate for an equatorial 



passage of the compact object is indeed predicted to be 
well above the Eddington limit. 

After the outburst changes in the emission lines persist 
for at least three years, with Fe II lines stronger than before 
and He I, Hen, and Nil lines weaker. The timescale for the 
extended Fell decay, at least, is similar to the expected 
viscous timescale of the disc of hundreds to thousands of 
days, so this may indicate the gradual recovery of the disc 
to its equilibrium state. As the system does not yet appear 
to have stabilised continued monitoring is important to 
determine if the system eventually recovers to the pre- 
outburst state or if it settles to a different level. 
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